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ABSTRACT 

We present line-strengths and kinematics from the central regions of 32 galaxies with 
Hubble types ranging from E to Sbc. Spectral indices, based on the Lick system, 
are measured in the optical and near infra-red (NIR). The 24 indices measured, in 
conjunction with models of the effects of varying abundance ratios, permit the breaking 
of age/metallicity degeneracy and allow estimation of enhancements in specific light 
elements (particularly C and Mg). The large range of Hubble types observed allows 
direct comparison of line-strengths in the centres of early-type galaxies (E and SO) with 
those in spiral bulges, free from systematic differences that have plagued comparisons 
of results from different studies. Our sample includes field and Virgo cluster galaxies. 
For early-type galaxies our data are consistent with previously reported trends of Mg2 
and Mgb with velocity dispersion. In spiral bulges we find trends in all indices with 
velocity dispersion. We estimate luminosity-weighted ages, metallicities and heavy 
element abundance ratios (enhancements) from optical indices. These show that bulges 
are less enhanced in light (a-capture) elements and have lower average age than early- 
type galaxies. Trends involving age and metallicity also differ sharply between early 
and late types. An anti-correlation exists between age and metallicity in early types, 
while, in bulges, metallicity is correlated with velocity dispersion. We consider the 
implications of these findings for models of the formation of these galaxies. We find 
that primordial collapse models of galaxy formation are ruled out by our observations, 
while several predictions of hierarchical clustering (merger) models are confirmed. 

Key words: galaxies: abundances - galaxies: evolution - galaxies: formation - galaxies: 
stellar content. 



C3 ' 1 INTRODUCTION 

The bulges of spirals, and the spheroidal components of ellip- 
tical (E) and SO galaxies exhibit many si milarities. These in- 



clude their colours and colour gradients ( Balcells & Peletier 



1994) as well as their morphological and kinematic proper- 
ties (Bender, Burstein & Faber 1993). The question then 
arises: do these similarities imply similar formation pro- 
cesses? Attempts to answer this question using morphology 
and kinematics have so far failed because the observed prop- 
erties are successfully reproduced by both models proposed 
for galaxy formation (i.e. primordial collapse, e.g. Carlberg 
1984 and hierarchical collapse, e.g. Kauffmann, White & 
Guiderdoni 1993 and Barnes & Hernquist 1996). On the 
other hand, the use of photometry to constrain the models 
is severely hampered by the degeneracy in colours with re- 
spect to age and metallicity (Worthey 1994, hereafter W94). 
Thus, in galaxy populations, the properties above place only 



* E-mail: rproctor@uclan.ac.uk; aesansom@uclan.ac.uk 



weak constraints on formation mechanisms and star forma- 
tion histories (SFHs). 

In order to provide a more sensitive tool for probing 
compo site populations, W 94 used the Lick/IDS spectral fea- 
tures (Faber et al. 1985), to estimate pseudo line-strengths 
(indices), for a range of single age, single metallicity stellar 
populations (SSPs). W94 demonstrated that, although in- 
dividual indices are affected by age/metallicity degeneracy, 
their sensitivities to both age and metallicity vary. Conse- 
quently, W94 was able to find combinations of indices that 
broke the degeneracy. Over the last few years, models of 
SSPs have been refined and new indices added. Vazdekis 
et al. (1996; hereafter V96) and Vazdekis (1999a,b; here- 
after V99) used more recent (and complete) isochrones for 
the SSP calculations, as well as including the Mgl and Ca 
triplet indices in the NIR (Diaz, Terlevich & Terlevich 1989). 
Higher order Balmer absorption lines (Rg and H 7 ) were also 
included. We shall refer to the collection of indices above as 
the Lick indices. 

Many studies of early-ty pe galaxies have measured 
small numbers of Lick indices (Worthey, Faber fc Gonzalez 
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1992; |Davies, Sadler fe Peletier 19931 |Gorgas et al. 1997 
Greggio 1997; [Tragcr ct al. 2000a ). These studies show 
that, when the commonly observed Mg2 or Mgb indices are 
plotted against indices centred on Fe features (e.g. Fe5227, 
Fe5335 or their average <Fe>), most early-type galaxies lie 
outside the range of values covered by SSPs based on so- 
lar metal abundance ratios. In these galaxies the Mg2 and 
Mgb index values are significantly larger than those of SSPs 
with similar Fe index values. This illustrates the well known 
Mg excess in early-type galaxies (e.g. Faber 1973; O'Connell 
1976; Worthey et al. 1992; Henry & Worthey 1999). Indices 
centred on C feat ures are also se e n to lie outside the rang e 
covered by SSPs ( |Worthey 1998| ; |Henry fe Worthey 1999TJ , 
implying C is another element whose relative abundance is 
enhanced in these populations. These studies also show that 
a positive correlation exists between central velocity disper- 
sion and certain Mg indices (e.g. Bender et al. 1993,1998 and 
Kuntschner 2000; hereafter K00). These trends have been 
likened to trends in colour-magnitude diagrams and inter- 
preted as a metallicity-mass relation (i.e. high metallicities 
are found in high velocity dispersion galaxies due to the abil- 
ity of deep potential wells to retain metals, e.g. Faber 1973; 
Kodama et al. 1998). However, this interpretation makes 
no allowance for either the age/metallicity degeneracy in 
these indices or their enhancements with respect to other 
indices. It is also important to note that other metallicity- 
sensit ive indices show weak or no trends with velocity disper- 
sion (Fisher Franx fe Illingworth. 1996; Trager et al. 1998 



j0rgensen 1999). 

Few studies have includ ed sp iral bulges. However, 
Jablonka, Martin fe Ari moto (1996) and Idiart, de Freitas 
Pacheco & Costa (1996) observed Lick indices in the cen- 



tres of bulges. They found correlations of the small number 
of indices observed with both bulge luminosit y and velocity 
dispersion. Goudfrooij, Gorgas & Jablonka ( |l99S| ) also re- 
port measurements of a small number of indices in bulges 
and suggest that the Mg and C excess observed in early 
t ypes is also present in bulges. However, in a previous study 
(Proctor, Sansom & Reid 200C), we found that a sample of 4 



spiral bulges lay closer to the solar abundance ratio SSPs, in 
the Mg2-<Fe> plane, than do early types. Unfortunately, 
none of the bulge studies a bove were emission corrected, 
and only Goudfrooij et al. (1999) was fully calibrated to 
the Lick system, making full comparisons with early-type 
galaxies difficult. This highlights one of the advantages of 
the present study in which such systematic uncertainties 
are significantly reduced by the identical instrumentation 
and reduction procedures used. 

While the sensitivities of Lick indices to abundances of 
individual elements, such as Mg and C, complicate their in- 
terpretation, they also provide a powerful tool for probing 
star formation histories in galaxies, using models of galac- 
tic chemical evolution (e.g. Vazdekis et al. 1997; Sansom & 
Proctor 1998). This is due to the difference between produc- 
tion sites of light elements and those of Fe-peak elements, i.e. 
while Mg is mainly produced in type II supernovae (SNII), 
Fe-peak elements are predominantly produced in type la su- 
pernovae (SNIa). Thus, if reasonable estimates can be made 
of Mg and Fe abundances, constraints can be placed on pos- 
sible SFHs. To quantify the effects of individual element 
abundance enhancements on Lick indices, Tripicco & Bell 
1995 (hereafter TB95) modelled the effects of doubling 10 



key elements in the synthetic spectra of 3 stellar types. They 
showed that Lick indices centred on C and Mg features (e.g. 
CNi, CN 2 or their average <CN>, C 2 4668, Mgi, Mg 2 and 
Mgb) are much more sensitive to the abundances of C and 
Mg than they are to overall metallicity ([Z/H]Q). These in- 
dices are then more sensitive to the abundance ratios of C 
and Mg (i.e. [C/Fe] and [Mg/Fe]) than to [Z/H]. On the 
other hand, the sensitivities to individual elements of the Fe 
indices (Fe4383, Fe4531, Fe5015, Fe5270, Fe5335 or Fe5406) 
and Ca indices (Ca4227 and Ca4455) are comparable to, 
or less than, their sensitivity to [Z/H]. It is this range of 
sensitivities to individual elements and [Z/H] that permits 
estimation of ages, metallicities and abundance ratios from 
index values, and thus the constraining of possible SFHs, in 
galaxy populations. 

In this paper we describe the measurement and analysis 
of 24 Lick indices in the centres of 32 galaxies with Hubble 
types ranging from E to Sbc. These galaxies were observed in 
a single observing run, with identical instrumental set-ups. 
We compare the indices and their correlations with veloc- 
ity dispersion in early and late-type galaxies. We use V99 
SSPs and the data of TB95 to estimate luminosity-weighted 
ages, metallicities and abundance ratios in our galaxy sam- 
ple. Finally the overall trends in our data are compared to 
the predictions of our galactic chemical evolution code for 
models of galaxy formation. 

In Section H the observations and data reductions are 
described, including calibrations and emission corrections. 
In Section ^ central values of indices and kinematics are pre- 
sented. Comparison of the trends in our data are made with 
previous observations. Luminosity-weighted age, metallicity 
and abundance ratio estimates are described in Section m 
In Section [| we discuss the interpretation of our results in 
terms of composite models of possible star formation his- 
tories. In Section ^| we give some discussion of our results 
and draw our conclusions. Appendix A gives details of ve- 
locity dispersion corrections. In future papers and in Proctor 
(PhD thesis - in preparation) the spatially resolved results 
will be reported and the data further exploited to recover 
more detail of the SFHs, using galactic chemical evolution 
models. 



2 OBSERVATIONS AND DATA REDUCTIONS 
2.1 Sample Selection 

The data presented here are from observations made dur- 
ing time awarded for two separate PPARC PATT propos- 
als. The first was to test the prediction that young ellip- 
tical galaxies should be devoid of hot gas (e.g. Ciotti et al. 
1991). Bright (B% < 13) early-type galaxies (—5 < T < -2), 
with X-ray emission temperatures and luminosities well con- 
strained by ROSAT observations, were selected. The sec- 
ond proposal detailed an investigation of the stellar popula- 



t We use the standard notation; [X/Y]=log[(X/Y)-log(X Q /Y Q )] 
where X and Y are the masses of either individual, or groups 
of, elements. The term abundance is used to refer to the mass 
of a specific element divided by that of H ([X/H]), while abun- 
dance ratios are specified with respect to Fe ([X/Fe]). Overall 
metallicity is represented by [Z/H]. 



© RAS, MNRAS 000, 



Stellar populations in galaxy spheroids. 3 



Table 1. WHT Observations. T type, Hubble type and half-light radius (r e ) are from de Vaucouleurs et al. 1991 (hereafter RC3). Central 
velocity dispersion (cto) and radial velocity (RV) are the values derived from our blue spectra and are for the central 3.6x1.25 arcsec. 
Estimated errors are given in brackets (see Section 2.5.1 for details of derived kinematics). Distances are mainly from Tully 1988 with 
the exception of more distant galaxies where radial velocities from RC3 were used. Exposure times and ISIS slit position angle (PA), 
which is normally along the minor axis of each galaxy, are given. Group membership is from Tully 1988. 
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* Not along minor axis (Section 2.2). 
t Distances calculated from radial velocities given in RC3 assuming Hq=75 km s — 1 Mpc -1 . 



tions along the minor axes of bright (By < 13) spiral bulges 
(0 < T < 4). Highly inclined galaxies (inclination > 75°) 
were selected to minimise the effects of disc contamination 
in the outer regions of the bulges. Highly inclined galaxies 
with prominent dust lanes covering the centre of the bulge 
were removed from the sample, as these would give little 
information about the bulge centres. Both studies aimed to 
use the same range of Lick indices to investigate the SFH of 
galaxies and to make estimates of luminosity- weighted ages, 
metallicities and abundance ratios at a number of points 
across the galaxies, to estimate gradients. Fortunately, the 
two observation runs were scheduled back-to-back allowing 
observations to be made with identical instrument set-ups. 
The two studies result in a sample of 32 galaxies with 24 
indices measured (Table |l|) . 



2.2 Observations 

Long-slit spectroscopic observations along the minor axes of 
11 Es, 6 SOs and 16 spiral bulges (bulges) were obtained be- 
tween 1998 February 28 and 1998 March 3 with the WHT on 
La Palma. The double beam spectroscope (ISIS) was used 
with a 5700 A dichroic and Tektronix 1024 square CCDs. 
On the blue arm a 300 line/mm grating was used giving 
a wavelength coverage of 3995-5495 A at a dispersion of 
1.5 A per pix el. Th is range covers 16 indices calibrated by 
Gorgas et al. ( 1993 ) including the extensively observed Mg2 



and <Fe> i ndices, as well as the recentl y calibrated H<5 and 
H7 indices (Worthey & Ottaviani 1997). On the red arm a 



600 line/mm grating was used, giving a wavelength cover- 
age of 8275 A to 9075 A at a dispersion of 0.8 A per pixel. 
This range covers the Mgl feature at 8807 A and the highly 
metallicity-sensitive Call triplet (Diaz et al. 1989; V96). The 
plate-scale was 0.36 arcsec per pixel on both red and blue 
arms. The slit, of length 4 arcmin and width 1.25 arcsec, 
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was placed along the minor axes of the galaxies with the 
exception of NGC 2831 and NGC 2832 which, being in close 
proximity on the sky, were observed simultaneously; i.e. with 
a position angle of 45° in both galaxies. A maximum expo- 
sure time of 1500 seconds was adopted to facilitate cosmic 
ray removal. Multiple exposures of individual galaxies were 
obtained to achieve the desired signal-to-noise (giving index 
errors of approximately 5% at r e /2). Seeing was ~1.5 arcsec. 
The total exposure time and position angle of the slit are 
given in Table [l]. Tungsten lamp exposures, for flat-fielding, 
were obtained each night on the blue arm. However, due to 
known fringing effects on the red arm, red tungsten lamp 
exposures were taken just before or after every object ex- 
posure, with the telescope tracking the object. A total of 5 
flux calibration standards and 24 stars (from Faber et al. 
1985) for calibration of the Lick indices were observed. A 
neutral density filter (ND1.8) was used in the stellar obser- 
vations. Observations of faint calibration stars and tungsten 
lamp exposures were obtained, with and without the neu- 
tral density filter, to allow removal of the filter's spatial and 
spectral responses. The sample of Lick calibration stars was 
selected to possess index values spanning the range of values 
expected in our galaxy sample. The calibration star sample 
was also chosen for good overlap with stars with H £ and H7 
measurements reported by Worthey & Ottaviani (1997]), as 



1997b, as 
lex (piaz 



well as stars used in the calibration of the Call index (I 
et al. 1989). All calibration stars possess known heliocentric 
radial velocities. 



2.3 Basic reductions 

Unless otherwise stated, data reductions were carried out 
using the CCDPACK, FIGARO and KAPPA packages of 
Starlink software. Bias removal was carried out by the sub- 
traction of an average bias frame, normalised to the average 
value in the over-scan region, in each object frame. After 
conversion from electrons to photons, variance arrays were 
generated and propagated throughout the reductions. In the 
blue, flat-fielding was achieved by division of target frames 
by the normalised average of tungsten lamp exposures ob- 
tained on the same night. However, on the last night suit- 
able tungsten lamp exposures were not obtained. For this 
night the flat-field from the first day was used (the day for 
which arc exposures were most similar). Division of target 
frames by tungsten lamp exposures leaves the spectra bi- 
ased by the smooth spectral response of the lamp. This is 
removed at flux calibration. However, during the flat-fielding 
procedure, features in the blue tungsten lamp spectra were 
identified that moved independently of wavelength calibra- 
tion. The features were in the range 4000 - 4600 A and were 
identified with features in the dichroic response. The effects 
of these features are included in our statistical errors as de- 



tailed in Section 2.4.2. Flat-fielding of the red data was car- 
ried out by division of each target frame by the normalised 
average of the bracketing tungsten lamp exposures. Stellar 
frames were divided by the normalised neutral density fil- 
ter response. Cosmic rays and bad rows were removed by 
interpolation across the affected areas. Wavelength calibra- 
tion was carried out by comparison with arc lamp exposures 
taken just before and/or after each exposure. An accuracy 
of better than 0.1 A was consistently achieved in both red 
and blue calibrations. All object frames were extinction cor- 



Table 2. The resolution of the Lick calibrations are given for 
reported indices. For all indices, with the exception of G4300, off- 
sets represent average differences between the published index val- 
ues and our observations for 24 Lick calibration stars. For G4300 
the offset was modelled as a function of the observed value (see 
text). N.B. <Fe> is the average of Fe5270 and Fe5335 while CaT 
is the sum of Ca2 and Ca3. 
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* Features affected by dichroic response. 

** Typical error given by Diaz et al. (1989) as 5%. 

tLick offset=4.34-0.749G4300ra m (see Section Ea2) 



rected using the extinction curve appropriate for La Palma. 
Flux calibrations derived from multiple observations of sin- 
gle stars varied by less than 1.5% across the region of CCD 
used, while those derived from differing stars varied by less 
than ~ 5%. All frames were flux calibrated using the aver- 
age of the calibration curves of 5 flux calibration stars. Sky 
estimates were made using the outermost regions of the slit 
that were not significantly vignetted. After sky subtraction 
galaxy frames were co-added to form a single frame for each 
galaxy. The spiral galaxy NGC 4100 was found to be dom- 
inated by emission and was excluded from further analysis. 
Due to the presence of telluric absorption lines above 8920 
A , reliable NIR indices could not be determined for galaxies 
with recession velocities above 2200 km s" 1 (i.e. NGC 2831, 
2832, 5353, 5354, 5908 and 5987). 

2.4 Calibrations using stellar spectra 

The original calibrations of stellar line-strengths with pho- 
tospheric parameters, upon which the SSPs used here are 
base d, were carried ou t using data from the Lick/IDS scan- 
ner ( Faber et al. 1985 ). This instrument has a spectral res- 
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olut ion that varies with wavelength (Worthey & Ottaviani 
1997|)7]rhe Lick data were also not flux calibrated. There- 
fore, in order to calibrate our index measurements to the 
Lick system, it is necessary to compensate for both the dif- 
ferences in flux calibration between our data and that in 
the calibration data, as well as the differences in spectral 
resolution achieved. 



2.4-1 Correction for spectral resolution 

For each index, the value of the Lick spectral resolution (o\l 
in Table Q) was estimated from Fig. 7 of Worthey & Ot- 



taviani (1997), at the wavelength of the mid-point of the 
central band. The spectra of our sample of Lick calibration 
stars were then broadened to the appropriate calibration res- 
olution (ol), for each index, by convolution with a Gaussian 
of width ob given by: 



2 

O-B 



2 



2 



(1) 



The instrumental broadening for our data (07) was es- 
timated from arc lines and found to be 1.5 A in the blue 
and 0.7 A in the NIR. After appropriate broadening, stel- 
lar indices were evaluated using our own code. Wavelength 
range definitions supplied by Worthey on his home page were 
used. For the NIR indices, band definitions an d cali bration 
resolution (ctl) were taken from Diaz et al. ( 1989 ) (after 
allowing for a typographical error). Our code was tested us- 
ing the stellar data also provided on Worthey's home page 
for this purpose. Differences between the values given by 
Worthey and those derived by our code from the provided 
spectra were < 0.03 A for the line features, and < 0.002 mag 
for molecular band indices. These discrepancies are smaller 
than differences caused by re-calibration of Worthey's data 
to our wavelength resolution and are probably the result of 
differences in the handling of partial bins. 



2-4-2 Flux calibration correction 

In order to compensate for the differences in flux calibra- 
tion between our data and the stellar calibration spectra, 
the difference between measured and published values was 
calculated, for each index, in each of the observed calibra- 
tion stars. For all indices except G4300, we found no signifi- 
cant correlation between these differences and the measured 
values. Therefore, for all indices except G4300, the average 
difference is used as a final correction to the velocity disper- 
sion corrected values. For G4300, the differences between 
measured (G4300 rat i,) and published values exhibited a cor- 
relation with G4300 r(lU j given by: 



Offset = 4.340 - 0.749 x G4300 ratJ 



(2) 



The final correction to the velocity dispersion corrected 
values of G4300 was, therefore, calculated using this equa- 
tion. However, for one galaxy (NGC3769), G4300,-a™ lay sig- 
nificantly outside the range of values covered by our stellar 
sample. Consequently, the value of G4300 in this galaxy, 
while reported here, was omitted from further analysis. 
Comparison of the measured indices in the calibration sam- 
ple with the published data is shown in Table |j| For blue 
indices with band definitions above 4600 A, the RMS scat- 
ter about the offset in our sample is generally dominated 



by the RMS error associated with individual Lick observa- 
tions. However, the scatter in our data is significantly greater 
than this error for indices with band definitions below 4600 
A (with the exception of Ca4227). This is the wavelength 
range affected by t he p oor removal of the dichroic response 
identified in Section 2.3. For indices in this wavelength range 



the excess scatter (calculated in quadrature) of our data 
compared to the Lick error has been included in our error 
calculations. This turns out to be a conservative error es- 
timate for all indices but Ca4227, whose error we assume 
is underestimated. For the NIR indices differences between 



the stellar data of Diaz et al. (1989) and our stellar data 
were used for flux calibration correction. We note that the 
scatter in our data is greater than t he RMS error per single 
observation in the Diaz et al. (1989) data (given as typically 
5%). However, these indices are not used for the purposes of 
absolute age/metallicity estimates and this uncertainty has 
not been included in our errors. For all indices, the error in 
calibration to the Lick system was calculated as the stan- 



dard error in the stellar data, i.e. 
number of calibration stars. 



where N=24 is the 



2.5 Analysis of galaxy spectra 

To derive accurate indices from a galaxy spectrum it is first 
necessary to obtain accurate estimates of recession velocity 
and velocity dispersion from the spectra. This allows the red- 
shift to be taken into account and the indices to be corrected 
for velocity dispersion (using the polynomials detailed in 
Appendix A). Gala xy da ta must also be corrected for flux 
calibration (Section 2.4.2) and emission. 



2. 5. 1 Measurement of galaxy kinematics 

Measurements of galaxy kinematics were carried out on both 
the red and blue data using the Fourier quotient technique 
within the IRAF software package. This technique was used 
as the associated statistical errors are significantly less than 
those associated with the cross-correlation technique at the 
low velocity dispersions typical of spiral bulges. For the 
purposes of velocity dispersion analysis, H/3 and [OIII]5007 
emission lines were removed by linear interpolation across 
affected regions in galaxies showing strong emission. Veloc- 
ity dispersion and recession velocity were estimated using 
each of the calibration stars. Final values and errors were 
taken as the average and standard error in the individual 
estimates. Typical errors for velocity dispersion and reces- 
sion velocity were found to be ~ 3.1 km s" 1 and ~ 2.6 km 
s _1 respectively. A comparison of our results for velocity dis- 
persion with the average values in Prugniel & Simien (1997) 
is shown in Fig. ^. In this figure the average value and RMS 
scatter in the individual measurements quoted in Prugniel & 
Simien (1997) are taken as the velocity dispersion and its er- 
ror respectively. Despite differences in spectral and spatial 
sampling, agreement is reasonably good between the data 
sets, with a one-to-one line having a reduced \ 2 of 1-7. Our 
galaxy indices were corrected for velocity dispersions as de- 
scribed in Appendix A. 
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sonable agreement (~ ± 0. 15 A ). The results of our OIII 



estimates are giv en in Table 3.1. 

Osterbrock shows that the line-strengths of H<5 



and H7 in emission are less than H/3, with line-intensities 
relative to H/3 of approximately 0.25 and 0.5 respectively 
across a large range of conditions. However, the continuum 
level in the spectra of all our galaxies show a reduction of 
approximately 50% between H/3 and H7. There fore, using 
our estimates of [OIII] emission, the Gonzalez (1993) cor- 



relation between [OIII] and H/3 emission, and Osterbrock 
data for H7 and H<5, we applied the following correc- 



tions: 



Fe5015 = Fe5015 ratl ,+AFe5015, 
H/3 = H/3 raw +AH/3, 
H7 = H 7row +AH7, 
RS = B6 raw +AHS, 



AFe5015 = -[OIII] 
AH/3 = -0.7[OIII] 
AH7 =-0.7[OIII] 
AH<5 = -0.35[OIII] 



Figure 1. Comparison of measured central velocity dispersions 
with average values from Prugniel & Simien (1997) for 22 galaxies 
common to both studies. The one-to-one line is shown. Prugniel 
&; Simien (1997) error was estimated as the RMS scatter in the 
individual measurements included in the database which include 
a large range of measurement methods and apertures. Errors in 
our data are the standard errors in values derived from individual 
calibration stars for each galaxy. 



2.5.2 Emission Correction 

The H/3, H7 and H<5 indices suffer from line-filling in galax- 
ies exhibiting emission. Fe5015 is also affected by [OIII]5007 
emission i n such galaxies, while Mgb is affe cted by [NIJ5199 
emission (Goudfrooij fc Emsellem 1996). We estimated 



[OIII] 5007 emission in all galaxy spectra in an effort to com- 
pensate for the effects of emission in H/3, 7, S and Fe5015. No 
attempt was made to estimate the (relatively small) correc- 
tions to Mgb. Follo wing a procedure similar to that used 
by Gonzalez (1993), each galaxy spectrum was divided by 
a series of template spectra. These templates were made by 
red-shifting and broadening each star used in the blue kine- 
matic analysis to the measured galaxy values. The aim of 
dividing the galaxy spectrum by a well matching template 
is to remove common spectral features around the [OIII] line 
prior to the measurement of its equivalent width. The band 
definitions used for our [OIII] line index are as follows: 

Side band 1: 4990.0-5001.0 A 
Centre band: 5001.0-5011.0 A 
Side band 2: 5011.0-5022.0 A 

The average of values derived using each stellar tem- 
plate was taken as the [OIII] index for each galaxy. The 
use of the average of a group of well matching stellar spec- 
tra differs from the method used by Gonzalez of creating 
individual templates for each galaxy. Also, [OIII] band defi- 
nitions differ between the two methods. Gonzalez estimated 
[OIII] emission in four galaxies common with our sample 
whil e thr ee of the galaxies were observed by Kuntschner et 
al. ( 2001 ). Differences in slit width and orientation make 
direct comparison difficult. However, the results are in rea- 



It can be seen that, due to the reduction in continuum 
level, H7 is as sensitive to emission in absolute terms (A) as 
H/3, while the H5 feature is only half as sensitive. However, 
due to the range of strengths of both H<5 and H7 in SSPs, 
estimates of age and metallicity made using these indices 
are significantly less affected by emission than those made 
using the H/3 index. 

An H/3 emission index similar to the [OIII] emission in- 
dex described above was also defined as a check on the [OIII] 
index. Comparison of [OIII] and H/3 emissi on in dices showed 
that most galaxies follow the Gonzalez (1993) correlation 
well. However, six galaxies were noted (all spiral bulges) 
with significantly aberrant behaviour. Three bulges (NGC 
4157, NGC 4217 and NGC 4312) show H/3 emission sub- 
stantially greater than that expected from the [OIII], while 
three (NGC 3254, NGC 3769 and NGC 4313) show strong 
[OIII], but with no detectable H/3 emission. These late-type 
galaxies do not follow the Gonzalez correlation. H/3 values 
for these six galaxies are omitted from our analysis. 



3 RESULTS 

3.1 Central Values 

Central values for velocity dispersion, given in Table [j], are 
those derived from the central 3.6x1.25 arcsec of the blue 
observations with centres defined as luminosity peaks. Veloc- 
ity dispersion values derived using blue spectra are greater 
than those derived from NIR spectra by an average of ~ 
10 km s _1 with scatter about this of ~ 20 km s _1 . How- 
ever, velocity dispersion profiles in some galaxies (e.g. NGC 
4192) differ significantly between blue and NIR data. This 
suggests that the two wavelength ranges may be detecting 
differing kinematic populations. 

Cen tral index values and errors for blue indices are given 
in Table O. Values for Fe5406 in NGC 2831 and NGC 2832 



were not determined due to the high recession velocity in 
these galaxies redshifting the blue side-band of this index 
outside the observed spectral range. Indices are corrected for 
both velocity dispersion and emission and converted to the 
Lick system. Reduction and calibration errors have all been 
included in the quoted errors. Also included, for indices with 
band definitions below 4600 A, are the er rors d ue to poor re- 
moval of the dichroic response (Section 2.4.2). Comparison 
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Table 4. NIR central index values (first line) and errors (sec- 
ond line). All reduction and calibration errors are included in the 
errors. Omitted galaxies are those whose high recession veloci- 
ties redshift one side-band of these indices into the region of the 
spectrum affected by telluric lines. 



Galaxy Cal 



Ca2 



Ca3 



Mgl 









(A) 




(A) 




(A) 




(A) 


NGC 


2549 


1 


774 


1 


528 


3 


878 





759 









078 





076 





062 





042 


NGC 


2683 


1 


697 


1 


201 


3 


685 





665 









101 





092 





079 





057 


NGC 


2832 


















NGC 


2831 












" 




" 


NGC 


3226 


1 


737 


4 


242 


3 


483 





910 









117 





110 





100 





071 


NGC 


3254 


1 


687 


1 


175 


3 


650 





671 









116 





106 





097 





068 


NGC 


3301 


1 


623 


1 


355 


3 


775 





575 









094 





088 





074 





051 


NGC 


3607 


1 


639 


1 


169 


3 


552 





806 









091 





089 





079 





059 


NGC 


3608 


1 


797 


4 


263 


3 


634 





854 









100 





096 





086 





060 


NGC 


3623 


1 


711 


1 


181 


3 


658 





748 









092 





086 





072 





051 


NGC 


3769 


1 


745 


3 


965 


3 


561 





550 









300 





268 





277 





222 


NGC 


4157 


1 


573 


1 


073 


3 


730 





543 









188 





167 





161 





131 


NGC 


4192 


1 


966 


1 


458 


1 


170 





652 









095 





089 





075 





051 


NGC 


4203 


1 


797 


1 


276 


3 


573 





864 









121 





111 





101 





069 


NGC 


4216 


1 


794 


1 


159 


3 


743 





732 









086 





082 





070 





046 


NGC 


4217 


2 


427 


1 


827 


3 


877 





770 









574 





518 





576 





446 


NGC 


4291 


1 


652 


4 


197 


3 


429 





679 









105 





104 





096 





075 


NGC 


4312 


1 


530 


3 


535 


3 


989 





344 









596 





518 





523 





331 


NGC 


4313 


1 


734 


1 


622 


4 


145 





710 









124 





112 





102 





072 


NGC 


4365 


1 


900 


1 


430 


3 


699 





723 









109 





105 





098 





074 


NGC 


4374 


1 


740 


1 


109 


3 


438 





597 









104 





102 





095 





075 


NGC 


4419 


1 


854 


1 


143 


3 


813 





630 









106 





099 





088 





060 


NGC 


4526 


1 


910 


1 


311 


3 


793 





823 









116 





108 





099 





068 


NGC 


4552 


1 


603 


3 


771 


3 


257 





705 









111 





105 





097 





074 


NGC 


4636 


1 


629 


4 


009 


3 


539 





707 









128 





119 





111 





081 


NGC 


4697 


1 


660 


1 


343 


3 


625 





703 









110 





102 





092 





062 


NGC 


5322 


1 


858 


4 


471 


3 


635 





778 









110 





105 





097 





072 


NGC 


5354 


















NGC 


5353 


















NGC 


5746 


1 


882 


4 


403 


3 


852 





812 









116 





107 





096 





067 


NGC 


5908 


















NGC 


5987 



















Figure available from 
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Figure 2. Comparison of our results for early-type galaxies to 
published data from Trager (1998) for commonly quoted indices. 
Also shown as triangles are results for Fe5270 and Fe5335 from 
two galaxies from Davies et al. (1993). 



of the results of our index measurements with galaxies in 
common with Trager (1998) are shown in Fig. |[ Two galax- 
ies from the study of Davies et al. (1993) are included for 
the Fe indices. The velocity dispersion sensitive Fe5270 and 
Fe5335 indices agree within errors. However, for Mg2, an off- 
set of ~ 0.008 mag is observed. Such an offset would cause 
metallicity and abundance ratio estimates, made using this 
index alone, to vary by ~ 0.04 dex. Part of the offset is prob- 
ably the result of the slightly wider (1.4 arcsec) and longer (4 
arcsec) aperture used in the Lick/IDS observations reported 
in Trager's thesis. Increasing the size of the central region 
from 3.6x 1.25 arcsec to 4.5x 1.25 arcsec (equal surface area) 
in our measurements results in an average reduction in Mg2 
of ~ 0.002 mag while Fe5270 and Fe5335 are both reduced 
by ~ 0.01 A. We therefore conclude that, for these indices 
at least, our data contain no significant systematic biases 
compared to previous observations. 

In Table [| values and errors for the fully calibrated NIR 
calcium triplet and Mgl features are presented. Values are 
from the central 3.6x1.25 arcsec based on luminosity peaks 
in the NIR. Quoted errors include reduction, calibration and 
velocity dispersion errors. Velocity dispersion correction was 
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Figure 3. Mg2 and Mgb with central velocity dispersion (ao in 
km s _1 ). Circles show our E and SO galaxies. The thick line shows 
the mean trend from K00. Thinner lines show trends from other 
authors as detailed in the text. N.B. Both Trager et al. (2000b) 
and K00 estimated correlations on a log(index) scale, resulting in 
the curvatures in these plots. 



based on velocity dispersion values derived from the NIR 
galaxy spectra. 



3.2 Correlations between kinematics and indices 

3.2.1 Early-type galaxies 

Correlations between indices sensitive to a-elements (e.g. 



served in early-type g 


alaxies (Bender et al. 1993; Gonzalez 


1992; j0rgensen 1997 


; |Bernardi et al. 1998; IConcannon et 


al. 2000; Kuntschner 2000; Trager et al. 2000b). Fig. |3|sum- 



marises results for Mg2 and Mgb from these studies. Al- 
though there is generally good agreement for the slopes of 
the correlations, there are offsets between studies. Such off- 
sets may result from differences in calibration to the Lick 
system, differences in aperture size and/or orientation and 
systematic differences in velocity dispersion estimates. Val- 
ues of indices in the early-type galaxies of our sample are 
also shown in Fig. ^. Correlations from the fully calibrated 
K00 study of Fornax cluster galaxies (shown as thick lines in 
Fig. |^), most closely match our results, although our sample 
covers a significantly narrower range of velocity dispersion 
than the K00 sample (a result of our selection of bright, 
nearby ellipticals). 

In Fig. ^ we show plots of selected metallicity-sensitive 
indices from the blue spectra against log velocity disper- 
sion for our galaxy sample. Indices omitted are those most 
severely aff ecte d by the dichroic response problem outlined 
in Sections 2.3 and [2,4.2| . Table |E| details the fits (minimis- 
ing x~) °f indices against velocity dispersion for both early 
and late-type galaxies in our sample. Values given are from 
fits with y-axis errors only (which dominate in most cases). 
Errors in slope and intercept are taken as half the difference 
between the fit using index errors for x~ minimisation and 
that using velocity dispersion errors. From Table ^ it can 



Table 5. Correlations between indices and velocity dispersion. 
The number of galaxies (N) in the fit, slopes, intercepts, (un- 
weighted) correlation coefficients (r) and x 2 for best fit correla- 
tions are given. Derivation of errors (given in brackets) are de- 
scribed in the text. For Mg2 and Mgb in early type galaxies, the 
correlations of K00 are also given. N.B. For Mgb the K00 cor- 
relation has been transformed from magnitudes (as presented in 
K00) to Angstrom used in this work. NGC4313 was omitted from 
correlations for late-type galaxies (see Section B.2.; 



Index 


N 


dI/d(log <T ) 


Intercept 


r 


x 2 


EARLY TYPE GALAXIES 








Ca4227 


17 


0.22 (9.36) 


0.76(22.03) 


0.21 


97 


Fe5015 


17 


-0.50(30.18) 


7.17(71.08) 


0.02 


105 


Fe5270 


17 


-0.84 (3.13) 


5.18 (7.38) 


-0.16 


108 


Fe5335 


17 


-0.62 (3.10) 


4.32 (7.30) 


-0.18 


66 


<Fe> 


17 


-0.72 (2.50) 


4.73 (5.90) 


-0.19 


137 


Fe5406 


15 


-0 38 CI 27) 


2.95 (3.00) 


-0.21 


18 


CNi 


17 


0.086(0.855) 


-0.089(2.013) 


0.19 


21 


CN 2 


17 


0.113(0.741) 


n i c\7(^ 7AK\ 


u. zo 


9fi 


<CN> 


17 


0.099(0.792) 


-0.098(1.866) 


0.21 


47 


C 2 4668 


17 


-1.99(13.45) 


12.73(31.68) 


-0.12 


218 


Mgi 


17 


0.066(0.216) 


0.003(0.510) 


0.35 


438 


Mg 2 


17 


0.114(0.178) 


0.062(0.419) 


0.51 


392 


K00 Mg 2 


13 


0.191(0.023) 


-0.127(0.054) 






Mgb 


17 


2.68 (1.76) 


-1.24 (5.39) 


0.61 


245 


K00 Mgb 


13 


2.60 (0.40) 


-1.06 (1.36) 






KS A 


17 


-2.11(20.71) 


2.61(48.70) 


-0.24 


29 


R5 F 


17 


-1.26 (3.32) 


3.05 (7.80) 


-0.43 


10 


H 7 A 


17 


-2.83(44.85) 


1.04(105.53) 


-0.19 


29 




17 


-1.66(31.03) 


2.66(73.04) 


-0.21 


10 


H/3 


17 


-1.74 (0.66) 


5.68 (1.56) 


-0.71 


19 


CaT 


13 


-2.35 (1.94) 


13.31 (4.51) 


-0.53 


52 


Mgi 


13 


-0.28 (0.86) 


1.40 (2.00) 


-0.51 


16 


LATE TYPE GALAXIES 








Ca4227 


14 


1.68 (0.22) 


-2.44 (0.48) 


0.90 


21 


Fe5015 


14 


4.51 (0.34) 


-3.77 (0.75) 


0.95 


21 


Fe5270 


14 


3.03 (0.14) 


-3.46 (0.30) 


0.93 


19 


Fe5335 


14 


2.86 (0.57) 


-3.50 (1.26) 


0.92 


35 


<Fe> 


14 


2.95 (0.24) 


-3.49 (0.52) 


0.96 


10 


Fe5406 


14 


2.09 (0.23) 


-2.65 (0.48) 


0.97 


24 


CNi 


14 


0.337(0.088) 


-0.687(0.190) 


0.92 


6 


CN 2 


14 


0.347(0.093) 


-0.669(0.202) 


0.91 


8 


<CN> 


14 


0.340(0.089) 


-0.683(0.193) 


0.91 


10 


C 2 4668 


14 


11.80 (1.53) 


-18.78 (3.39) 


0.96 


84 


Mgi 


14 


0.230(0.042) 


-0.377(0.091) 


0.94 


155 


Mg 2 


14 


0.413(0.027) 


-0.623(0.060) 


0.97 


170 


Mgb 


14 


5.37 (0.22) 


-7.43 (0.48) 


0.98 


70 


HJ A 


11 


-12.34 (1.93) 


25.72 (4.18) 


-0.93 


24 


H<5 F 


14 


-6.04 (0.08) 


13.80 (0.17) 


-0.95 


38 




14 


-14.33 (2.00) 


26.41 (4.40) 


-0.88 


24 


F~7f 


14 


-7.04 (1.63) 


14.68 (3.59) 


-0.84 


10 


H/3 


9 


-1.44 (0.74) 


5.09 (1.64) 


-0.81 


16 


CaT 


12 


-1.03 (5.88) 


10.25(12.06) 


0.28 


50 


Mgi 


12 


0.39 (0.31) 


-0.12 (0.63) 


0.78 


9 
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Figure 4. Central metallicity-sensitive indices plotted against log of central velocity dispersion (ctq in km s -1 ). Thick lines are fits to our 
spiral bulge data (solid symbols with symbol size largest for SO/a and smallest for Sbc). Open symbols show early- type galaxies (circles 
for Es and squares for SOs). Thin lines show trends from K00 for early-type galaxies. Error bars include all reduction and calibration 
errors. Uncertainties in calibration to the Lick system are shown as an isolated error bars in bottom right of each plot. 



be seen that in the metallicity sensitive indices of early-type 
galaxy sample we detect no significant (3a) slopes. How- 
ever, our results for Mg2 and Mgb are consistent with the 
positive trends with velocity dispersion noted by previous 
studies (Fig. ^|). If we compare the scatter of our Mg2 and 
Mgb data about our best fit lines (0.02 mag and 0.3 A re- 
spectively) with their scatter about the K00 correlations 
(Fig. ^), we find no significant difference. We tested the fit 
of the K00 correlations to our early-type galaxy data for 
all common indices. For the majority of indices our scat- 
ter about the K00 correlations is similar to scatter of the 
K00 data. Consequently, in Fig. ^| we use K00 correlations 
(shown as thin lines) for all common indices for purposes of 
comparison to bulges. Most previous authors report weak 
or no trends in Fe indices (Fisher el al. 1996; Trager et al. 
1998; J0rgensen 1999), with velocity dispersion for early- 
type galaxies. Within the narrow velocity dispersion range 
covered by our sample, we echo this finding (Fig. ^| and Table 
However, K00 found correlations for all the metallicity- 
sensitive indices he observed. We note the possibility that 
such correlations would have been present in our data had 
the velocity dispersion range of our study been larger. For 
Fe5406 we find values consistently above the K00 correla- 
tion. However, the side bands of this index lie in a region of 



our spectra where vignetting effects are beginning to appear. 
Consequently, it is possible that a small undetected system- 
atic effect may be biasing this index in our results. Previous 
authors have reported that the calcium indices (Ca4227 and 
Ca4455) follow similar trends to Fe ind ices, rather than the 
trends in enhanced indices such as Mg2 ( Vazdekis et al. 1997 ; 
Trager et al. 199c). We again echo this finding within the 



narrow range of velocity dispersions of our sample. 

Plots of age sensitive indices with velocity dispersion 
are shown in Fig. [| In early-type galaxies no significant 
trends with velocity dispersion are evident in the H<5 and 
H7 indices. However, for the H/3 index a ~3<j correlation is 
found (Fig. |B| and Table |B|) . The slope of this weak correlation 
is steeper than that found by K00 (dH/3/d(log cr)~ -0.4) 
which is shown as a thin line in Fig. |E| The presence of such 
a correlation in the highly age sensitive H/3 index, while the 
more degenerate H<5 and H7 indices show no correlations, 
suggests the presence of a trend of increasing age and/or 
decreasing metallicity with velocity dispersion in the early- 
type galaxies. 

SO galaxies, shown as open squares in Fig. ^, follow sim- 
ilar trends to elliptical galaxies in all indices. The similarity 
in indices in E and SO types can also be seen in Fig. ^ where 
the values of key indices are compared for all Hubble types. 
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Figure 6. Central indices with Hubble (T) type. Open symbols 
are early types and filled symbols are late types. 



In this plot indices show no correlation with Hubble type 
for spiral bulges. 



Figure 5. Age sensitive indices against log of central velocity 
dispersion. Symbols and errors as Fig. 4. 



3.2.2 Late-type galaxies 

All metallicity-sensitive indices in the blue spectra of our 
late-type galaxy sample show strong correlations with veloc- 
ity dispersion (Fig. Table [s]), including those not shown 
in Fig. ^. One outlier to these trends for late-type galaxies is 
NGC 4313, which exhibits a strong central depression in log 
a. This galaxy has been omitted from the line fitting proce- 
dure. The highest velocity dispersion bulges show index val- 
ues coincident with those of early-type galaxies. Indeed, for 
all indices in our sample, early and late types form a contin- 
uous locus in the index-velocity dispersion plane. However, 
the slopes of the correlations in late-type galaxies are sig- 
nificantly steeper than those found by K00 for early-type 
galaxies (Fig. In general, the scatter about the correla- 
tions of late-type galaxies is also smaller than that of early 
types. Anti-correlations with velocity dispersion are also evi- 
dent among the late-type galaxies in all age sensitive indices 
(Fig. |^) with low velocity dispersion bulges having stronger 
hydrogen absorp tion l ines (aberrant emission galaxies de- 
tailed in Section 2.5.2 are omitted from line fitting in the 



case of H/3). The strong trends of both metallicity and age 
sensitive indices with velocity dispersion in late-type galax- 
ies suggests that SFH is closely associated with the depth of 
potential well in the centres of bulges. 
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Table 6. Index values from correlations at log(<ro)=2.0. The per- 
centage difference between the values from the two correlations is 
also given. 



Figure available from 
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Index 


K00 


Spirals 


% 


C 2 4668 


5.37 


4.81 


-10% 


Fe5015 


5.02 


5.25 


+5% 


Mgi 


0.114 


0.083 


-27% 


Mg 2 


0.255 


0.194 


-24% 


Mgb 


4.14 


3.31 


-20% 


Fe5406 


1.69 


1.52 


-10% 


<Fe> 


2.69 


2.40 


-10% 



Figure 7. CaT and Mgi against central velocity dispersion. Sym- 
bols and errors as Fig. 4 



3.2.3 Red data 

CaT (the sum of Call at 8498 A and 8542 A) and Mgi are 
plotted against velocity dispersion in Fig. pj Trends with 
velocity dispersion are not as clear in these indices. In early- 
type galaxies, the Mgi index shows behaviour with veloc- 
ity dispersion different to that of the Mgi, Mg2 and Mgb 
indices in the blue. However, it is difficult to make direct 



comparisons because, as previously suggested (Section 3.1), 
NIR indices may be sampling different stellar populations 
to those sampled by blue indices. 



3.2.4 Potential disc contamination of late-type galaxies 

Late and early-type galaxies in our sample exhibit similar 
index values in the region where their velocity dispersions 
overlap (150 - 200 km s _1 ). However, towards lower disper- 
sions, bulges lie systematically further from the K00 trends 
for early- type galaxies (Fig. W). It was a concern that this 
behaviour may be the result of di sc con tamination. Khos- 
roshahi, Wadadekar & Kembhavi (2000) carried out bulge- 
disc decomposition of a number of edge-on (i>50°) spiral 



galaxies. They show that for Hubble types earlier than Sbc 
(T<4) the bulge-to-disc central luminosity ratio always ex- 
ceeds 10; i.e. no more than 10% of the observed light origi- 
nates in the disc. They also show that this ratio increases to 
~ 1000 in the case of Sa galaxies. The worst-case assumption 
for metallicity-sensitive indices is that of a smooth contin- 



uum contribution from the disc. Under such circumstances, 
a maximum reduction of ~ 10% is expected for both line 
and molecular band indices. We compare the K00 correla- 
tions for metallicity-sensitive indices in early-type galaxies 
with our correlations for late types in Table H. Values of in- 
dices on the correlations are given, for our bulge data and the 
K00 correlations, at log ao=2.0 (mid-range value for bulges). 
Percentage differences between the two correlations can be 
seen to vary from -27% to +5%. Clearly the difference can 
not be modelled by the simple addition of continuum to the 
bulge light. It should also be noted that, for the Mg indices, 
a minimum bulge contribution of ~ 25% is required to ac- 
count for the differences between correlations. This is well 
in excess of the 10% maximum expected from Khosroshahi 



et al.(200C), and the difference increases at velocity disper- 
sions less than 100 km sec - . We also note that there is no 
evidence that the indices of late-type galaxies correlate with 
Hubble type (Fig. ^|) as might be expected if significant disc 
contamination were present. 

We therefore conclude that while metallicity-sensitive 
indices could be depressed by as much as 10% by disc con- 
tamination, this effect is not evident in our data, nor can 
it explain the observed differences between correlations in 
early and late- type galaxies. 



3.3 Diagnostic index plots 

In this section we compare the results of our index deter- 
minations to the values predicted for SSPs. We use V99 
SSPs for our analys is. Th ese (web published) SSPs are based 
on Bertelli et al. ( 1994) isochrones (as detailed in V96), 
and Vazdekis (1999a) transformations to the observational 
plane. We use V99 SSPs rather than those of W94 and 



Worthey & Ottaviani (1997), as V99 make use of th e mor e 



up-to-date and complete isochrones of Bertelli et al. ( 1994 ) . 
V99 also includes the CaT and Mgi indices in their calcula- 
tions of SSP index values. Metallicities derived from degen- 
eracy breaking diagnostic plots using V99 SSPs are gener- 
ally higher by ~ 0.1 dex than those implied by W94 SSPs. 
For galaxies older than ~ 5 Gyr, V99 SSPs also imply ages 
younger by ~ 0.15 dex. 

The <Fe> vs Mg2 diagnostic plot is shown in Fig. ^| 
This plot clearly shows the enhancement of the Mg2 index 
(or, equivalently, the deficiency in the <Fe> index) in early- 
type galaxies with respect to solar abundance ratio SSPs. 
The plot also demonstrates the age/metallicity degeneracy 
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Figure 8. <Fc> vs Mg2 (left) for our WHT data and three galaxies from Proctor et al. 2000 (solid squares). Values for fourth galaxy 
from Proctor et al. 2000 (NGC 3623) agree (within errors) with values from the WHT observation. All other symbols are as Fig. 4. Grid 
lines show V99 solar abundance ratio SSPs at ages of 1.5, 2, 3, 5, 8, 12 and 17 Gyr and [Fe/H] SS p of -1.7, -1.0, -0.5, -0.25, 0.0, +0.25, 
and +0.5 dex. In the left plot the dashed line shows the correlation of Gorgas et al. (1997) for early-type galaxies, which is consis tent 
with our data. The thick solid line shows the range of index values covered by the primordial collap se m odels described in Section 5.1, 
while the short, thinner solid line shows the range of values covered by the merger models in Section 5.3, Error bars in the bottom right 
of each plot refer to errors incurred in calibrating onto the Lick system. 



of these two indices. The grid includes SSPs ranging from 1.5 
to 17 Gyr in age and with metallicities ([Fe/H]ssp) cover- 
ing the range [Fe/H]ssp=-L7 to +0.5. The good agre emen t 
between our data and the correlation of Gorgas et al. (1997) 
for early-type galaxies is also illustrated. Late-type galaxies 
seem to span the SSP grid. This echos the tentative finding 
of our previous study (Proctor et al. 2000), that late-type 
galaxies exhibit abundance ratios closer to solar than those 
found in ell iptical galaxies. The 4 bulges observed in Proctor 
et al. (200C) were not fully calibrated to the Lick system due 
to a lack of observations of suitable calibration stars. How- 
ever, we note that most fully calibrated studies find that 
small positive corrections are required to compensate Mg2 
for the lack of flux calibration in the Lick observations, while 
Fe5270 and Fe5335 require only very small corrections (e.g. 
Worthey & Ottaviani 1997; K00). The Lick offsets given in 
Table H for the pres ent data set were thus applied to the 
Proctor et al. ( |2000| ) data. The Sa galaxy NGC 3623 was 
observed in both studies. The two independent values for 
both Mg2 and <Fe> in this galaxy were within errors. Con- 
sequently, this galaxy was not included when the results of 
the Palomar study were added to Fig. | (solid squares). Ad- 
ditional errors of 0.01 mag and .04 A have been added in 
quadrature to Proctor et al. (200C) errors for Mg2 and <Fe> 
respectively to allow for uncertainty in the corrections to the 
Lick system. 

Given the tendency, reported by previous authors, for 
Ca indic es to follow similar trends to Fe indices ( Vazdekis et 
al. 199 /|; (Trager et al. 199g| ), the NIR CaT and Mgl indices 
allow plotting of a diagnostic plot similar to <Fe> — Mg2 in 
the blue. This diagnostic plot (CaT— Mgl) is also shown in 
Fig. ^. In this figure, early-type galaxies lie within with the 
SSP grid. A simplistic interpretation would be that this in- 
dicates [Mg/Ca]=0 and, assuming that Ca follows the same 
trends as Fe, a value [Mg/Fe]=0. However, when these in- 
dices are plotted against age sensitive indices (e.g. H/3) prob- 



lems with the Mgl index become apparent (Section 3.3.1) 
and this interpretation can not be sustained. 

3.3.1 Breaking the degeneracy 

Plots that most clearly break the age/metallicity degeneracy 
are shown in Figs ^| and |l(j. These plots show metallicity- 
sensitive indices against H/3 and compare galaxy values to 
V99 SSP predictions. We show indices against H/3 since it is 
the most age sensitive index and other age sensitive indices 
(R5a,f, G4300 and H7a,f) are among those effected by the 
dichroic response problem in our data. 

The H<5 and H7 indices were also not among the indices 
whose sensitivities to element abundance ratios were mod- 
elled by TB95. Fig. [] shows indices sensitive to Fe on the 
left, while indices sensitive to abundance ratios are shown 
on the right. Early- type galaxies in Fig. |^ exhibit different 
trends in Fe and abundance-ratio-sensitive indices. In these 
galaxies, the increasing H/3 with line strength in Fe indices 
contrasts with the decreasing H/3 with line strength in abun- 
dance ratio sensitive indices. The trend in Fe indices suggests 
that Fe abundance is anti-correlated with age in early-type 
galaxies, while the growing disparity between Fe and abun- 
dance ratio sensitive indices suggests increasing abundance 
enhancement in light elements with age. 

In late-type galaxies, H/3 decreases with all metallicity- 
sensitive indices in the blue (Fig. ^). The position of these 
galaxies with respect to the SSP grids is also more consistent 
between plots of Fe indices and those sensitive to abundance 
ratios. This reflects abundance ratios closer to solar in these 
objects, as noted in Proctor et al. 2000. The positions of 
late-type galaxies also suggests that the populations are rel- 
atively young (luminosity-weighted ages < 5 Gyr) with a 
wide range of metallicities. 

Fig. hd shows a plot with H/3 of the blue Ca4227 (top 
plot). Comparison of this plot with Fig. [] shows the simi- 
larity in behaviour of Ca and Fe indices previously reported 
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Figure 9. Mctallicity-sensitive indices plotted against H/3. Symbols as Fig 4. Fe indices are shown on the left, while indices sensitive to 
Mg and C abundance ratios are shown on the right. Uncertainty in calibration to the Lick system is shown top right of each plot. Short, 
thick and thinner lines show primordial collapse and merger models respectively. These and the grid lines are as in Fig. pi i-i i — i 
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wavelength bands are reasonably small and uniform across 
our galaxy sample, we attempt to interpret these diagrams. 

The CaT index again shows the similarity in behaviour 
of Ca and Fe indices. That is, the trend of increasing H/3 
with CaT, is similar to those exhibited by blue Fe indices. 
However, in the MgI-H/3 plot, while some galaxies retain the 
behaviour of blue Mg indices, many galaxies are displaced to 
extremely low metallicities when compared to the SSP grids. 
These galaxies tend to be both the oldest, and the highest 
velocity dispersion. However, Mgl is a weak index, prone 
to measurement problems, therefore further observations of 
Mgl are needed to ascertain its behaviour relative to SSPs. 
We therefore draw no strong conclusions from either the 
MgI-H/3 or Mgl-CaT plots. 



4 ESTIMATION OF AGE, METALLICITY AND 
ABUNDANCE RATIO ENHANCEMENT 

Our aim is to use the indices of solar-neighbourhood abun- 
dance ratio SSPs as the basis for estimating luminosity- 
weighted ages, metallicities and abundance ratios of galax- 
ies. This requires us to estimate the effects on indices of 
the non-solar abundance ratios observed in many galaxies 
(e.g. Davies et al. 1993). Thus we aim to construct grids 
of SSPs with varying age, [Fe/HJssp and abundance ra- 
tio for comparison with observations. This requires knowl- 
edge/estimation of: 



• The abundance ratio pattern in the local stars used to 
construct SSPs. 



• The difference between the local abundance ratio (SSP) 
pattern and that in the galaxy populations being studied. 

• The effects such a difference would have on photospheric 
conditions (surface gravity (log g), effective temperature 
(T e //) and luminosity) in whole populations i.e. the effects 
on isochrones. 



• The effects of differences in abundance ratios on the 
strength of individual indices in stars, assuming fixed pho- 
tospheric conditions. 



Figure 10. Calcium sensitive indices; Ca4227 (top) and NIR CaT 
(middle) plotted against H/3. Also shown is the Mgl index (bot- 
tom) against H/3. Uncertainty in calibration to the Lick system is 
shown bottom right of each plot. Symbols as Fig. 4. 



flVazdekis et al. 1997| ; [Trager et al. 1998| ) with Ca4227 sug- 
gesting lower metallicities at increased age. Also shown in 
Fig. nfj are the NIR CaT and Mgl indices plotted against 
H/3. A note of caution here is that we are comparing in- 
dices derived from NIR spectra with H/3 derived from the 
blue spectra. As previously noted, we suspect that the two 
wavelength ranges may be sampling differing populations. 
However, under the assumption that the differences in age 
and metallicity of the populations sampled by NIR and blue 



Armed with the above it is possible to estimate cor- 
rections to the indices of solar abundance ratio SSPs for 
non-solar abundance ratios, for comparison to observations. 
The above points are discussed in turn below. 



4.1 The local abundance pattern 

Stellar metallicity estimates, used in the construction of 
both W94 and V99 SSPs, were based on analysis of solar- 
neighbourhood stars (e.g. Edvardsson et al. 1993 in V99 
and Hansen & Kjaergaard 1971; Gustafsson, Kjaergaard & 
Andersen 1974 in W94). The wavelength ranges used gen- 
erally contain a large number of Fe lines, thus provid- 
ing good estimates of iron abundance ([Fe/H]). However, 
solar-neighbourhood stars possess considerable scatter in the 
abundance ratios of individual elements. In high metallic- 
ity stars ([Fe/H] > solar) abundance ratios show moderate 
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Table 7. Mass of individual elements as a fraction of total met- 
als (Xj/Z)© in the Sun. Data is from Cox (2000). The total mass 
fraction of metals in the Sun (Zq) is assumed to be 0.0189. Ele- 
ments enhanced in low-metallicity, solar-neighbourhood stars and 
in galaxy populations are identified by a + . These are the elements 
which are included in the enhanced group (E) in each case. The 
final column indicates (with Y) the elements modelled by TB95. 
The last row in the second column shows the total fraction of 
all tabulated elements in the Sun. For the other 3 columns totals 
indicate solar proportions of elements identified by + or Y. 



A[Z/H] = A[Fe/H] + AA[E/Fe] 



(4) 



Element 


(Xi/Z) s 


Low [Fe/H] 


Galaxy 


Modelled 






stars 


populations 


by TB95 


C 


0.1619 




+ 


Y 


N 


0.0583 




+ 


Y 





0.5054 


+ 


+ 


Y 


Ne 


0.0921 


+ 


+ 




Na 


0.0018 


+ 


+ 


Y 


Mg 


0.0343 


+ 


+ 


Y 


Al 


0.0030 


+ 


+ 




Si 


0.0370 


+ 


+ 


Y 


S 


0.0193 


+ 


+ 




Ar 


0.0054 


+ 


+ 




Ca 


0.0034 


+ 




Y 


Cr 


0.0009 






Y 


Fe 


0.0719 






Y 


Ni 


0.0039 








Total 


0.9986 


0.7017 


0.9185 


0.8749 



scatter about solar valu es ( Edvardsson et al. 1993 



Feltz- 



ing & Gustafsson. 1998). Thus for SSPs, the calibration of 



which averages a large number of local stars, we may assume 
[Z/H] SS p=[Fe/H]ssp for [Fe/H]ssp>0. On the other hand, 
solar- neighbourhood stars with [Fe/H] < exhibit non-solar 
abundance patterns, with a-elements (e.g. O, Ne, Na, Mg, 
Al, Si, S, Ar, Ca) enhanced with respect to Fe p e ak elements 
(e.g. Fe, Ni , Cr) flRyan Norris 



Bessell 1991 



et al. l! )93 



Feltzing & Gustafsson. 1998; tdiart & Thevenin 



Edvardsson 



200C). The a-element abundances are also enhanced with 
respect to C, which has solar abundance ratio (or below) 
down to very low metallicities (Ryan et al. 1991). Abun- 
dance ratios of a-elements ([a/Fe]) in these stars increase 
linearly from [a/Fe]~0.0 at [Fe/H]=0.0 to [a/Fe]~+0.3- 
0.5 (depending on the element concerned) at [Fe/H]= — 1. 
Moderate scatter in the abundance ratios of individual ele- 
ments about these trends is again observed. As SSPs with 
[Fe/H] ssp < are therefore based on stars with non-solar 
abundance ratios, we assume that [Fe/H] ssp more closely 
approximates [Fe/H] in the stars than [Z/H]. Thus, in low 
metallicity SSPs (and any other non-solar abundance ratio 
population) we must assume a relationship between [Z/H], 
[Fe/H] and [E/Fe] (where E is the mass of all elements with 
enhanced abundance rat ios - see Table colu mn 3). As non- 
linear effects are small (Tantalo, et al. 1998), and following 



Trager et al. (2000a; hereafter TOOa), we assume a linear 
relationship: 



[Z/H] = [Fe/H} + A[E/Fe]. 
The differential form of which: 



(3) 



may be used to derive estimates of the factor A in low metal- 
licity stars by comparing a solar composition with one in 
which all a-elements are doubled in abundance ([Fe/H]=0, 
[E/Fc]=+0.301). From the values given in Table M it can 
be seen that a-elements constitute 70% of metals at solar 
abundance ratio. Doubling all these elements therefore gives 
A[Z/H]=+0.25 (assuming AY=2.2AZ; Lebreton et al. 1999 
and an internally consistent handling of H). This yields an 
estimate of A — 0.83. Non-linear effects and uncertainties in 
the handling of He result in an uncertainty in the value of 
A of approximately ± 0.06. For solar neighbourhood stars 
we characterise the trend in [E/Fe] with [Fe/H] in the range 
-l<[Fe/H]<0 by: 

[E/Fe] = -0.45[Fe/H}. (5) 

The value 0.45 in Equation |^ is the overall enhancement in 
a-element abun dance s with respect to solar in the data of 
Salasnich et al (200C) for low metallicity field stars. Com- 
bining Equations ^ and |B| we get: 

[Z/H] =0.63[Fe/H]. (6) 

We use this equation to estimate [Z/H]ssp in SSPs with 
-l<[Fe/H]ssp <0. 



4.2 Abundance ratio patterns in galaxies 

Given the well known over-abundance of Mg with respect 
to Fe in elliptical galaxies (e.g. O'Connell 1976; Worthey 
et al. 1992) one approach to modelling abundance ratios 
in galaxies would be to use the pattern observed in low- 
metallicity, solar-neighbourhood stars which show a similar 
enhancement in Mg. However, while galaxy studies show 
Mg sensitive indices (Mgi, Mg2 and Mgb) to be enhanced 
with respect to Fe indices such as Fe5227, Fe5335 (Gorgas 
et al 1997; Vazdekis et al. 1997; K00; Trager et al. 2000a; 
Fig. ^| this paper) , many of these studies show indices cen- 
tred on Ca features (Ca4227, Ca4455) to be un-enhanced 
(see also Fig. [l(]). This is not the only difference between 
local, low-metallicity stars and galaxy populations. C sen- 
sitive indices (CNi, CN2 and C24668) are also enhanced 
in galaxy populations (Vazdekis et al. 1997; K00; Fig. ^ 
this paper). The differences between abundance patterns in 
low-metallicity, solar-neighbourhood stars and that in high- 
metallicity galaxy populations are interesting, but not sur- 
prising given the difference in both metallicity and environ- 
ment. They do, however, render this approach to modelling 
galaxy abundance ratios unworkable, as many Lick indices 
are particularly sensitive to C abundance. Here we make the 
assumption that enhancement in Mg abundance reflects an 
equal enhancement in the abundances of all a-elements in 
Table (?], with the exception of Ca which is assumed to fol- 
low Fe peak element abundances. C is also assumed to be 
enhanced. For galaxy population abundances we have there- 
fore defined two groups; the 'Fe-like' elements (Ca, Cr, Fe 
and Ni) and the 'enhanced' elements (C, N, O, Ne, Na, Mg, 
Al, Si, S and Ar - see Table ^, column 4) with the abundance 
of each element, in each group, enhanced by the same factor. 
This is similar to model 4 of T OOa, where a similar analysis 
was carried out (see Section 4.5). The two groups of elements 
combined represent 99.86% of the mass of metals present in 
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the solar photosphere (Table (5]) . Using these assumptions we 
hope to obtain reasonable estimates of age, metallicity and 
the degree of enhancement in galaxy populations. 

For the 3 galaxies with estimated [Fe/H]<0 it is nec- 
essary to compensate for the non-solar abundance ratios in 
the stars used to calibrate the SSPs. We use Equation |H| to 
calculate the stellar enhancement. This is added to the es- 
timated enhancement (relative to SSP). Finally, Equation ^ 
and the new value of [E/Fe] are used to obtain [Z/H]. 

4.3 Non-solar abundance ratio isochrones 

Until recently, understanding of the effects of non-solar 
abundance ratios on stellar populations was poor. How- 
ever, desire for a better description of the evolutionary 
tracks of globular clusters led to the development of theoret- 
ical isochrones for a-element enhanced populations. Salaris 
& Weiss (1998) and VandenBerg et al. (2000) modelled 
isochrones for the abundance pattern in metal-poor field 
stars which possess e nhan cements in a-element abundances 
They show that, at low metallici- 



as detailed in Section 4.1 



ties ([Z/H]<0), a-enhanced isochrones are well estimated by 
solar abundance ratio isochrones of the same [Z/H]. How- 
ever, they also show that at the highest metallicities in their 
studies ([Z/H]~solar) this no longer holds. Salasnich et al. 
(2000) modelled a-enhanced isochrones at metallicities simi- 
lar to those observed in galaxy centres ([Z/H]~ 0—0.5). They 
show that, at these metallicities, due to low atomic numbers 
and high ionisation potentials, a-elements make relatively 
small contributions to opacity. Consequently, high metallic- 
ity, a-enhanced isochrones possess a higher T e / / than solar 
abundance ratio isochrones of the same [Z /H] . The resultant 
isochrones are therefore best approximated by solar abun- 
dan ce rat io isochrones of lower [Fe/H] . Fig. 10 of Salasnich et 
al. ( 2000 ) shows that an isochrone for a 10 Gyr old popula- 
tion, with solar metallicity and an a-enhancement of +0.45 
(i.e. [Fe/H] = —0.37, Equation is best modelled by a so- 
lar abundance ratio isochrone of metallicity [Fe/H]= —0.25. 
At [Z/H]=+0.6 ([Fe/H]=+0.23) the effect becomes even 
more pronounced such that the a-enhanced isochrone strad- 
dles the solar abundance ratio isochrones of [Fe/H]=+0.0 
and [Fe/H]=+0. 3. Co nsequently, an approximate use of the 
Salasnich et al. (200C) isochrones is to assume that the ad- 
dition of a-elements to high metallicity populations results 
in no change in the isochrones. Equivalently, the position of 
the isochrones may be considered to depend upon [Fe/H] in 
the population rather than [Z/H]. This is cl early only an ap- 
proximation. However, the Salasnich et al. (2000) isochrones 



suggest that, at high metallicities, it is better than assuming 
a-enhanced isochrones to be best modelled by solar abun- 
dance ratio with the same [Z/H]. Differences in abundance 
pattern between metal-poor field stars and galaxy popula- 
tions are not expected to make a significant difference to 
the behavior of the isochrones as the majority of significant 
contributors to opacity are amongst the elements assumed 
to show similar abundance ratios in both populations (O, 
Ne, Na, Mg, Al, Si, S, Ar and Ni), while no other individ- 
ual elements, including C and Ca, contribute significantly to 
opacity. 

As SSP indices are based on empirical relations for line- 
strengths versus T r jj, lo g g and [Fe/H] in stars, the results 
of Salasnich et al. ( |2000|) have a significant impact on the 



interpretation of SSPs. If it is assumed that isochrone posi- 
tions (e.g. T e //, log g) are dependent on [Z/H], then SSPs 
with [Fe/H] SS p=[Z/H]p p (where [Z/H] PO p is [Z/H] in the 
galaxy) provide the best match to enhanced populations. If, 
on the other hand, it is assumed that isochrone position is 
dependent on [Fe/H], then [Fe/H]ssp=[Fe/H]pop provides 
the best match. 



4.4 Estimating indices in non-solar abundance 
ratio stars 

TB95 modelled the effects of individually doubling the abun- 
dances of 10 key elements in the synthetic spectra of 3 key 
stellar types (cool giant, turn-off star and cool dwarf), as 
well as doubling all elements simultaneously. All models were 
evaluated at fixed T e // and log g, with values based on a 5 
Gyr, solar metallicity isochrone. TB95 assumed no change in 
the opacity distribution function when doubling individual 
elements. This approximation is appropriate as individual 
elements contribute little to the opacity. For the doubling of 
all elements, TB95 again assumed fixed T e // and log g, while 
an opacity distribution function appropriate for twice solar 
metallicity was used. The results are presented as variations 
of the 21 Lick indices modelled, in terms of a standard er- 
ror, caused by doubling each of the 10 elements. The results 
of doubling all elements are presented in the same manner. 
Indices modelled by TB95 include 16 present in our study. 
TB95 did not, however, model US, H7 Mgl or CaT. One of 
the elements modelled by TB95 was Ti. The abundance of 
this element is enhanced in low metallicity stars while it's 
atomic mass lies in the range of the Fe-like elements. Con- 
sequently, the decision as to whether to include Ti among 
the enhanced or Fe-like indices is difficult. We have therefore 
not included this (low abundance) element in our analysis. 
TB95 is the only published study of this type to date. 

The effects of changing abundance ratios can be esti- 
mated from the TB95 data by defining Ri,x', the fractional 
change in the ith index in the TB95 arrays (e,g, CNi, CN2, 
Ca4227 etc) when the abundance of the element X (e.g. C/H, 
N/H) is doubled. The value of the enhanced index (/;) can 
then be estimated from the solar abundance ratio value (F) 
by (following TOOa): 



92 



(7) 



Where Ex is the change in abundance of element 
X, i.e. Ex=A[X/H]. Doubling all elements (column 14 of 
TB95; Tables 4 to 6) can be handled in the same way i.e. 
Ez=A[Z/H]. It should be noted that the doubling of all 
elements in TB95 did not model a simple doubling of the 
stellar metallicity as the TB95 calculations were made as- 
suming fixed photospheric conditions (T e ff, log g). 



4.5 Constructing non-solar abundance ratio SSPs 

It is possible to estimate the effects of non-solar abundance 
ratios in populations by modelling SSPs as combinations 
of the cool giant/turn-off/cool dwarf stellar types of TB95. 
The luminosity-weighted sum of the factors in Tables 4 to 
6 of TB95 allows estimation of Ri,x values for SSPs. These 
can then be combined, by use of Equation [?], to estimate 
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Table 8. Fractional changes (Ri,x) in index values of 53/44/3 mix when abundances of individual elements (e.g. C/H) are doubled. 
These data are derived from Tables 4 to 6 of TB95 and arc shown here to illustrate the sensitivities to individual elements. 



Index 


C 


N 





Mg 


Fe 


Ca 


Na 


Si 


Cr 


Ti 


z 


CNi 


1.75 


0.54 


-0.46 


-0.13 


-0.03 


-0.07 


-0.03 


0.11 


-0.10 


0.03 


0.40 


CN 2 


1.09 


0.35 


-0.29 


-0.10 


-0.03 


-0.05 


-0.02 


0.12 


-0.05 


0.03 


0.29 


Ca4227 


-0.32 


-0.05 


0.10 


0.00 


0.05 


0.31 


-0.01 


0.00 


-0.01 


0.00 


0.24 


G4300 


0.27 


0.00 


-0.06 


-0.02 


-0.04 


0.01 


-0.01 


-0.01 


-0.02 


0.05 


0.04 


Fe4383 


0.08 


-0.00 


-0.02 


-0.05 


0.20 


-0.03 


-0.01 


-0.04 


0.00 


0.02 


0.12 


Ca4455 


-0.05 


-0.01 


0.01 


-0.01 


-0.08 


0.00 


-0.01 


-0.00 


0.07 


0.03 


0.16 


Fe4531 


0.00 


0.01 


0.01 


-0.01 


0.03 


0.00 


0.01 


-0.05 


0.04 


0.11 


0.14 


C 2 4668 


1.90 


0.00 


-0.27 


-0.06 


0.03 


-0.01 


-0.01 


-0.12 


-0.02 


0.04 


0.36 


H/3 


0.03 


0.00 


-0.00 


-0.04 


0.00 


0.00 


0.01 


0.01 


-0.04 


0.00 


-0.00 


Fe5015 


-0.00 


0.00 


-0.00 


-0.10 


0.09 


0.01 


0.01 


-0.03 


-0.02 


0.08 


0.14 


Mgi 


0.78 


-0.00 


-0.11 


0.26 


-0.10 


-0.01 


-0.01 


-0.05 


-0.01 


0.02 


0.21 


Mg 2 


0.11 


-0.01 


-0.03 


0.23 


-0.04 


-0.01 


-0.01 


-0.03 


-0.01 


0.02 


0.14 


Mgb 


-0.17 


-0.01 


-0.00 


0.37 


-0.07 


-0.00 


-0.01 


-0.04 


-0.10 


-0.01 


0.08 


Fe5270 


0.07 


0.01 


-0.00 


-0.05 


0.11 


0.02 


-0.01 


-0.01 


0.01 


0.02 


0.14 


Fe5335 


-0.05 


-0.01 


0.01 


-0.04 


0.20 


0.00 


-0.01 


-0.00 


0.03 


0.02 


0.14 


Fe5406 


0.03 


0.01 


-0.00 


-0.01 


0.17 


-0.01 


-0.01 


-0.00 


0.05 


0.02 


0.15 



the sensitivities of indices in SSPs to enhancement of el- 
eme nts in the cho sen abundance ratio pattern. Trager et 
al. (2000a; 2000b) used a 53/44/3 percentage luminosity- 
weighted combination of the three stellar types to simulate 
SSPs. Although not detailed in their paper, the 53/44/3 
combination is consistent with values given by W94 for the 
relative contributions of these stellar types to 3-17 Gyr SSPs 
of [Fe/H]g5p ~0. The similarity in the sensitivities (in per- 
centage terms) of cool giant and turn-off stars make the re- 
sults fairly robust to any reasonable combination. We have 
therefore adopted the TOOa luminosity weightings through- 
out this paper for ease of comparison. Fractional changes 
in indices for doubling individual element abundances, and 
doubling [Z/H], for this combination of stellar types, are 
given in Table ^. The sensitivities of indices to individual 
elements divides them into 2 groups; those with sensitivities 
to individual elements significantly larger than their sensi- 
tivity to Z (CNi, CN 2 , G4300, C 2 4668, Mgi and Mgb) and 
those whose sensitivities to individual elements are similar 
to or less than their sensitivity to Z (the rest). It is inter- 
esting to note that the strong sensitivities of the first group 
are mainly to C and Mg, both of which are enhanced el- 
ements in galaxy populations (e.g. see Fig. ^j, right side). 
Consequently, in galaxies, these indices are dominated by 
abundance ratio effects rather than [Z/H]. 

The next step in constructing non-solar abundance ratio 
SSPs is to select the SSP whose isochrone best matches that 
expected in the population under study. It is to the indices 
of this SSP that the TB95 corrections for non-solar abun- 
dance ratios are applied. On the basis of the isochrone mod- 
els available at the time (salaris & Weiss 1998; VandenBerg 



et al. 2 )00), TOOa assumed that the isochrone position was 



governed by [Z/H], thus an SSP with [Fe/H] SS p=[Z/H] P op 
was selected. However, if w e acc ept the implications of the 
subsequent Salasnich et al. (200C) study - that isochrone po- 
sitions are dependent on [Fc/H] rather than [Z/H] - then we 
must select an SSP with [Fe/H]ssp = [Fe/H]p'OP- Given the 
uncertainties in this aspect of the modelling we tested three 



methods for applying the TB95 data to galaxy observations; 
the method used by TOOa, which assumes isochrone shape 
to be governed by [Z/H], and t wo me thods designed to be 
consistent with Salasnich et al. (2000) isochrones. 



4.5.1 TOOa method 

Given the assumption that isochrone positions are governed 
by [Z/H], TOOa pointed out that there is only one way to 
achieve the required element abundance ratio enhancements. 
This involves the reduction in the abundance of Fe-like ele- 
ments, while enhanced element abundances are (marginally) 
increased to maintain [Z/H]. Abundances of elements not 
modelled by TB95 are assumed to remain constant. TOOa 
used Equation ^|to derive: 



A\Fe/H] = —AA[E/Fe], A[Z/H] = 



(8) 



Where E now refers to all elements enhanced with re- 
spect to Fe-like elements in galaxy populations. The data 
in Table |^, column 4 leads to a value of A = 0.942 for 
the galaxy abundance ratio pattern using the TOOa method. 
From Equation || this method requires a value of Ex (in 
Equation [?]) given by: 

E x = (1- A)A[E/Fe], X = C, N,0, Mg, Na, Si (9) 
While for Fe-like indices: 



E x = -AA[E/Fe], X = Fe, Ca, Cr 



(10) 



As this method directly estimates age, [Z/H] and [E/Fe], 
Equation |^ is used to calculate [Fe/H]. The difficulty with 
this method is that it is based on the (now apparently false) 
assumption that, at the high metallicities observed in galaxy 
centres, isochrone positions are governed by [Z/H]. Nonethe- 
less, we have applied the TOOa analysis to our data using V99 
SSPs. 
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4-5.2 Methods based on Salasnich isochrones 



In light of the Salasnich et al. (2000) isochrones we chose to 
test two alternative approaches to the application of TB95 
data. As these methods are based on the assumption that, 
at the high metallicities present in the majority of galaxies, 
isochrone shape is governed by [Fe/H], we select the SSP 
with [Fe/H]ssp=[Fe/H]pop- We then ensure that the ap- 
plied enhancements involve no change in the abundance of 
Fe-like elements. We identify two ways of achieving these 
requirements with the TB95 data. The first is simply to in- 
crease the abundance of each of the elements thought to be 
enhanced in galaxy populations, while keeping other element 
abundances constant; the E+ method. For this method, the 
values of Ex, used in Equation [?], are given by: 



E X = A[E/Fe], X = C, N,0, Mg, Na, Si 



(11) 



For other elements and Z in Table ^ Ex=0. As this method 
is specifically designed to reflect high metallicity isochrones, 
it is not used for comparisons with observations of three low 
metallicity ([Z/H]<0) bulges. Instead, for these galaxies, we 
use the TOOa method which is consistent with the low metal- 
licity isochrones of Salaris & Weiss (1998) and VandenBerg 
et al. (2000). Metallicities and abundance ratios in these 
three g alax ies are transformed to solar scale as described in 



4.2 



Section 

The second approach is to double all elements (Z; Ta- 
ble |^) then reduce the abundances of Fe-like elements (Fe, 
Ca and Cr); the Fe— method. This method effectively dou- 
bles all elements except those in the Fe-like group. It should 
again be noted that doubling Z in the tables of TB95 does 
not represent a simple doubling of the metallicity in a real 
population, as the TB95 calculations were carried out at 
fixed T e f f and log g, i.e. with no movement of the isochrone. 
However, this is exactly the requirement of this method, as 
it assumes that addition of elements whose abundances are 
enhanced in galaxy populations leaves isochrone positions 
unchanged. For the Fe— method, the values of Ex used in 
Equation ^ are given by: 



E X = A[E/Fe], X = Z 
and 

E x = -A[E/Fe], X = Fe,Ca,Cr 



(12) 



(13) 



Again we have used the TOOa method for three bulges 
with [Z/H]<0. Despite the differences between TOOa and 
Fe— methods, there are strong similarities, as the fractional 
changes calculated for both methods are dominated by the 
reduction in Fe-like elements. However, the assumptions re- 
garding isochrones do result in significant diff erenc es be- 
tween estimated ages, as we shall see in Section 4.5.3. 



As both E+ and Fe— methods are based on SSPs of 
known [Fe/H], rather than [Z/H], we must re- write Equation 
| as: 

A[Z/H] = AA[E/Fe], A [Fe/H] = (14) 

The data in Table §, column 4 yields A = 0.941 for both 
E+ and Fe— methods. As both these methods directly esti- 
mate age, [Fe/H] and [E/Fe], Equation ^ is used to calculate 
[Z/H]. 
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Table 9. Comparison of Fe— , E+ and TOOa methods for mod- 
elling non-solar abundance ratio SSPs. Mean deviations from the 
best fit enhanced SSP models are given (see text for details). 



Index 



N 



Fe- 



E+ 



TOOa 







A/Toan 
iviL-elll 


v 2 
A 


Moan 

1V11_>CL11 


v 2 
X 


Moan 

lvlt_.Cl.il 


v 2 
X 


H<5 A 


32 


1.1 


74 


1.9 


160 


2.1 


179 


H<5 F 


32 


1.0 


81 


2.2 


228 


2.1 


212 


CNi 


32 


0.5 


32 


0.1 


25 


0.1 


21 


CN 2 


32 


0.7 


53 


0.5 


38 


0.1 


38 


Ca4227 


32 


-2.8 


378 


-3.3 


485 


-1.7 


194 


G4300 


32 


0.2 


52 


0.1 


49 


0.1 


51 


H 7 A 


32 


0.6 


49 


1.1 


79 


1.0 


80 


H7F 


32 


0.3 


22 


0.6 


32 


0,5 


29 


Fe4383 


32 


0.6 


26 


0.1 


19 


0,5 


21 


Ca4455 


32 


-1.0 


50 


-0.6 


22 


-1.0 


44 


Fe4531 


32 


0.2 


16 


0.3 


18 


0,1 


20 




39 
oz 


-u.o 


ou 


ft Q 

-u.y 


7-1 


ft 1 
-u. 1 


■ }J 


H/3 


26 


-1.4 


111 


0.5 


128 


-1.2 


139 


Fe5015 


32 


-0.6 


139 


1.2 


127 


-0.1 


109 


Mgi 


32 


-0.1 


102 


-2.1 


212 


-0.1 


97 


Mg 2 


32 


1.2 


128 


1.8 


187 


1.4 


152 


Mgb 


32 


-0.7 


89 


6.2 


1792 


-0.8 


98 


Fe5270 


32 


0,1 


67 


-1.4 


198 


0.7 


72 


Fe5335 


32 


0.2 


55 


-1.6 


163 


-0.5 


70 


Fe5406 


30 


1.9 


165 


-0.3 


24 


1.3 


91 


CaTt 


26 


-3.2 


(507) 


-2.8 


(376) 


-5.5 


(1189) 


MgF 


26 


6.8 


(1374) 


6,5 


(1287) 


6,1 


(1229) 


Total 






1769 




4060 




1756 



t Not included in x 2 minimisations. 



4-5.3 Comparison of results from different methods 

For all three methods, grids of non-solar abundance ra- 
tio SSPs were generated for [E/Fe] ranging from —0.3 
to +0.6 in steps of 0.025 dex. Grid spacings of 0.0125 
dex were used for 0.175<log(age)<1.225 and 0.025 dex for 
-0.5 <[Fe/H]< 0.75. Linear extrapolation of the V99 data 
from [Fc/H]=+0.4 to +0.75 was necessary due to the high Fe 
index values of a single SO galaxy (NGC 2549). Indices not 
modelled by TB95 (H<5 A , H5 F , R-y A and H 7f , Mgi and CaT) 
were assumed to have no sensitivity to abundance ratio. The 
best fit (minimum \ 2 ) was found for each galaxy, for each 
method. For the 6 galaxies observed to have aberrant emis- 
sion (Section 2.5.2), H/3 was omitted from the minimisation 
procedure. The Mgi and CaT indices were also omitted. 

Table ^ shows the average deviation (difference between 
observed and best fit non-solar abundance ratio SSP indices, 
as a multiple of our observational error) and \ 2 values for 
each index, over all observed galaxies, using each of the three 
methods. Comparison of best fits for E+ and Fe— methods 
shows the Fe— method to have a significantly lower total 
X 2 - We therefore dismiss the E+ method as too unrealistic. 
The TOOa method has marginally lower total \ 2 than the 
Fe— method. However, this difference hinges on a si ngle in - 
dex (Ca4227) whose error is underestimated (Section |2.4,2 ). 
If this index is excluded the Fe— method possesses a total 
X 2 value ~ 12% lower than the TOOa method. Comparison 



of the results for estimates of luminosity- weighted log(Age) , 
[Fe/H], [Z/H] and [E/Fe] from TOOa and Fe- methods are 
shown in Fig. [H]. It can be seen that for [Fe/H], [Z/H] and 



20 R.N. Proctor & A.E. Sansom 



Figure available from 
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Figure 11. Comparison of results obtained using the TOOa 
method to those using the Fe— method (see Section 4.5). Open 
circles are elliptical galaxies, open squares are SOs and filled sym- 
bols are late-type galaxies. 



[E/Fe] the results for the two methods are in good agree- 
ment. However, the TOOa method gives values of log(Age) 
significantly lower (by < 0.25 dex) than those derived by 
the Fe— method. This is a direct result of the difference in 
assumptions about the effects of non-solar abundance ratios 
on isochrone positions detailed in Section 4.5. The log(A ge) 
ordering agrees fairly well, to within ~ 0.1 dex. Table [10j 
presents results derived by the Fe— method for all galaxies 
wit h [Z/H ]>0, to be consistent with the recent Salasnich et 
al. ( |200C| ) isochrones. For galaxies with [Z/H]<0 the TOOa 
method has been applied in line with the low metallicity 
isochrones of Salaris & Weiss (1998) and VandenBerg et al. 
(2000). 

As a check on our decision to include all Balmer line 
indices (H/3, YL^a,f and U.8a,f) in the derivation of the val- 
ues given in Table [To], estimates were also made with dif- 
fering combinations of these 5 indices excluded from the fit- 
ting procedure. Table |ll] details the comparisons of these 
age/metallicity estimates with those given in Table [H]. Also 
shown are the scatters of the derived values about key cor- 
relations identified in Section |5| Offset and scatter in the 
derived values, when combinations of Balmer lines are ex- 
cluded from the fitting procedure, are small, with values 
<0.05 dex in cases when H/3 is not omitted. Indeed, even 
when all Balmer line indices are omitted, agreement with 
the values in Table hfl is relatively good. The comparisons 
presented in Table p3T therefore, show that our results are 
robust to the choice of Balmer lines to include in the fitting 
procedure. The results even suggest that, with a large num- 
ber of indices, reasonable estimates can be made without 
the Balmer lines. As many studies have been based on com- 
binations of indices such as H/3, Mg2 and <Fe>, we have 
also compared results from this combination in Table |l~ij 
We find only modest offsets. However, the scatter is large 
(~ 0.15 dex). Scatter about the key correlations are also 
found to be increased in this case. These results emphasise 



the advantage of the large number of indices included in this 
study. 

The effects on the results of ignoring enhancements are 
also shown in Table [HJ. To obtain this comparison, the best 
fits of our data (all 20 indices) to un-enhanced SSP index 
values were found. It was assumed that [Fe/H] ssp of the 
best fit SSP best represents the value of [Z/H] for the popu- 
lation. Comparison of the values of log(Age) and [Z/H], esti- 
mated in this way, to the values given in Table [Hj (derived by 
the Fe— method), show that the average difference between 
values are <0.05 dex in both parameters. Scatter about 
these differences are similarly <0.05 dex. We also tested 
the method employed by some authors (e.g. Gonzalez 1993; 
Kuntschner & Davies 1998) of combining the Mgb and <Fe> 
indices by taking their geometric mean {^/Mgb < Fe >) to 
make a new index ([MgFe]). Age and metallicity estimates 
are then made using this index and an age sensitive index; 
normally H/3. Results of the comparison between log(Age) 
and [Z/H] values derived using this method, and those de- 
rived from the Fe— method, are also shown in Table [n]. The 
average differences between the two sets of values is similar 
in magnitude to those found above for all 20 indices (with- 
out enhancement). However, the scatter in the differences is 
significantly greater (~ 0.15 dex). This result again empha- 
sises the benefit of including a large number of indices in the 
fitting procedure. 



4.6 Results from blue indices 

Values for luminosity- weighted log(Age), [Fe/H], [Z/H] and 
[E/Fe] derived by the Fe— method detailed above, are given 
in Table [u]. Plots of [Fe/HL [Z/H] and [E/Fe] against 
log(Age) are shown in Fig. Il2l To avoid confusion points 
with errors >0.2 dex, or limits, (see Table px| ) have been 
omitted from these plots. Confidence contours (1 a, allow- 
ing for 3 interesting parameters) are plotted on the right of 
Fig. [l2] as quadrilaterals with vertices at extremes projected 
onto the plane presented. The alignment of the contours in 
plots of [Z/H] and [Fe/H] against log(Age) clearly show the 
age/metallicity degeneracy. However, the size of the contours 
suggests that, for the majority of galaxies, the degeneracy 
has been broken. 

It can be seen in Fig. [li] that the ages of spiral bulges, 
SOs and Es form a continuous, overlapping sequence of in- 
creasing luminosity-weighted age, with bulges ranging from 
1.5 to 6 Gyr old, while SOs and Es range from 2 to 7 Gyr 
and 4 to 13 Gyr respectively. SOs also appear more Fe rich 
than Es. This is reflected in the strong central Fe features 
seen in SOs (e.g. Figs ^ and ^). 

We fitted lines (by x 2 minimisation) to both the early 
and late-type data plotted in Fig. The four galaxies iden- 
tified by asterisks in Table [jij were omitted from the fit- 
ting procedure. Table [l^ give the results of our fits. These 
have been calculated as the average of the fits obtained 
when y-axis and x-axis deviations are minimised separately. 
Errors quoted in Table |l2| are half the difference between 
the values from these two fits. Significant correlations are 
shown as lines in Fig. [l^. We find a strong anti-correlation 
between [Fe/H] and luminosity- weighted log(Age) in early- 
type galaxies. However, it was a concern that this may be 
heavily influenced by the low age, SO galaxy (NGC 2549) . We 
therefore re-calculated the correlation omitting this galaxy. 
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Figure 12. Luminosity-weighted [Fc/H], [E/Fe] and [Z/H] arc plotted against log(Age). Recall E represents the enhanced element group 
(Table 7, column 4). Data identified with an asterisk in Table 10 are omitted from this plot. Correlations for early-type galaxies are 
shown as lines. Key: »=l*u[g«s, □=S0, Q = E, x=Virgo cluster galaxies. 1 a confidence contours from our fits are shown on the right. 
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Table 10. Central values of luminosity- weighted log(Age), [Fe/H], [E/Fe] and [Z/H] for the 32 galaxies in our sample. Extent of 1 a 
confidence contours are given as errors (in brackets). Data points with an asterisk are omitted from Fig. 12. 



Galaxy 


Log(Age) 


[Fe/H ] 


[E/Fe] 


[Z/H] 


Ellipticals 


NGC2832 


0.875(0.069) 


0.175(0.050) 


0.300(0.025) 


0.457(0.074) 


NGC2831 


0.613(0.069) 


0.075(0.050) 


0.250(0.038) 


0.310(0.085) 


NGC3226 


1.025(0.056) 


-0.075(0.038) 


0.425(0.025) 


0.325(0.061) 


NGC3608 


0.950(0.044) 


0.125(0.025) 


0.275(0.013) 


0.384(0.037) 


NGC4291 


1.013(0.050) 


0.000(0.038) 


0.300(0.025) 


0.282(0.061) 


NGC4365 


0.988(0.044) 


0.175(0.025) 


0.250(0.013) 


0.410(0.037) 


NGC4374 


1.138(0.056) 


0.000(0.038) 


0.300(0.025) 


0.282(0.061) 


NGC4552 


0.988(0.044) 


0.150(0.038) 


0.325(0.025) 


0.456(0.061) 


NGC4636 


0.913(0.063) 


0.175(0.038) 


0.275(0.025) 


0.434(0.061) 


NGC4697 


0.712(0.075) 


0.300(0.025) 


0.200(0.025) 


0.488(0.049) 


NGC5322 


0.625(0.056) 


0.250(0.025) 


0.150(0.025) 


0.391(0.049) 


SOs 


NGC2549 


0.313(0.113) 


0.600(0.075) 


0.150(0.013) 


0.741(0.087) 


NGC3607 


0.750(0.056) 


0.250(0.025) 


0.250(0.025) 


0.485(0.049) 


NGC4203 


0.813(0.063) 


0.225(0.038) 


0.300(0.025) 


0.507(0.061) 


NGC4526 


0.587(0.038) 


0.325(0.038) 


0.175(0.025) 


0.490(0.061) 


NGC5354 


0.688(0.113) 


0.225(0.062) 


0.250(0.038) 


0.460(0.098) 


NGC5353 


0.863(0.081) 


0.300(0.050) 


0.275(0.025) 


0.559(0.074) 


Bulges 


NGC2683 


0.675(0.088) 


-0.025(0.050) 


0.200(0.038) 


0.163(0.085) 


NGC3254 


0.587(0.056) 


0.025(0.050) 


0.150(0.025) 


0.166(0.074) 


NGC3301 


0.338(0.050) 


0.175(0.050) 


0.075(0.013) 


0.246(0.062) 


NGC3623 


0.725(0.069) 


0.300(0.037) 


0.125(0.013) 


0.418(0.049) 


NGC3769 


<0.175(0.013)* 


-0.576(0.013)* 


0.234(0.113)* 


-0.356(0.346)* 


NGC4157 


0.512(0.169) 


-0.168(0.138) 


0.201(0.087) 


0.021(0.220) 


NGC4192 


0.300(0.056) 


0.350(0.038) 


0.150(0.025) 


0.491(0.061) 


NGC4216 


0.600(0.025) 


0.400(0.025) 


0.200(0.025) 


0.588(0.049) 


NGC4217 


0.463(0.363)* 


-0.050(0.438)* 


0.022(0.275)* 


-0.029(0.697)* 


NGC4312 


0.463(0.213)* 


-0.449(0.213)* 


0.227(0.237)* 


-0.235(0.436)* 


NGC4313 


<0.175(0.050)* 


0.325(0.050)* 


0.075(0.038)* 


0.396(0.085)* 


NGC4419 


0.313(0.044) 


0.050(0.050) 


0.175(0.038) 


0.215(0.085) 


NGC5746 


0.750(0.125) 


0.250(0.050) 


0.175(0.025) 


0.415(0.074) 


NGC5908 


0.475(0.069) 


0.350(0.075) 


0.125(0.025) 


0.468(0.099) 


NGC5987 


0.500(0.056) 


0.350(0.075) 


0.150(0.025) 


0.491(0.099) 



A slope of 0.672±0.193 was found indicating that this point 
is not having an excessive influence (cf Table We also 
detect a strong correlation between [E/Fe] and log(Age) in 
early- type galaxies. These tren ds are reflected in the index- 
index plots of Fig. ^ (Section 3.3.l[ ). The anti-correlation 
between [Z/H] (calculated by Equation]^) and log(Age) is a 
natural consequence of the somewhat stronger correlations 
with [Fe/H] and [E/Fe]. Both the (relatively young) SOs and 
the 5 Virgo cluster, ellipticals (which are amongst the old- 
est in our sample) appear to follow the same trends as the 
early-type sample as a whole. We find that elliptical galax- 
ies tend to be older and more metal-poor than SO galaxies. 
No significant correlations were found in late-type galaxies. 
However, the [E/Fe]— log( Age) plot shows that the data for 
early and late types form a continuous monotonic locus. This 
correlation may then be common to both Hubble types. 

Luminosity- weighted log(Age), [Fe/H], [Z/H] and 
[E/Fe] are plotted against log of central velocity dispersion 
in Fig. [13l Fits to these data (minimising \ 2 ) are given in 



Table [13[ Values given are from fits with y-axis errors only 
(which dominate in most cases). However, errors given are 
half the difference between the quoted fit and x-axis error 
only fit. Correlations with > 3<r significance are shown as 
lines in Fig. |l3[ There is a trend for luminosity-weighted 
age to increase with central velocity dispersion in our early- 
type galaxy sample. This trend is also consistent with the 
late-type data. However, we find no significance in the cor- 
relation of age versus velocity dispersion in late types (Table 
|l3|). The trend of increasing [E/Fe] with velocity dispersion 
for early-type galaxies suggested by Fig. jl3| also has very 
low significance particularly when compared to the strength 
of the correlation of [E/Fe] with log(Age) in early-types. 
Our data therefore suggest that, in the early- type galaxies 
in our sample, the stronger correlations of both [E/Fe] and 
velocity dispersion with log(Age) are the cause of the appar- 
ent trend of [E/Fe] with velocity dispersion. Bulges, on the 
other hand, while showing no significant correlations with 
age, show strong correlations of [Fe/H] and [Z/H] with veloc- 
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Table 11. Comparison of results of age/metallicity estimates omitting differing groupings of indices. Average offset and RMS scatter 
about offset (in brackets) are given for each of the derived parameters. Also shown (where appropriate) are the scatters about the most 
significant correlations detected (see Tables [l2| and |l3|) . N.B. 6 galaxies with aberrant emission (Section 2.5.2) were excluded from this 
analysis. 



Early Type Early Type Late Type 
Indices log(Age) log(Age) log(cr) 

omitted N Alog(Age) A[Fc/H] A[E/Fe] A [Z/H] vs [Fe/H] vs [E/Fe] vs [Fe/H] 



WITH ENHANCEMENT 



None 


20 


0.000(0.000) 


0.000(0.000) 


0.000(0.000) 


0.000(0.000) 


0.089 


0.041 


0.087 


H<5p,H7f 


18 


-0.020(0.026) 


0.005(0.016) 


0.008(0.021) 


0.012(0.020) 


0.089 


0.037 


0.077 




16 


-0.030(0.044) 


0.002(0.023) 


0.013(0.020) 


0.015(0.022) 


0.084 


0.031 


0.073 


H/3 


19 


0.021(0.066) 


-0.011(0.037) 


-0.004(0.026) 


-0.015(0.050) 


0.087 


0.050 


0.098 


H/3,H<5 a , f ,H 7j4 , f 


15 


-0.007(0.090) 


-0.009(0.042) 


0.006(0.027) 


-0.003(0.056) 


0.091 


0.059 


0.102 


All except 


3 


-0.070(0.171) 


0.004(0.129) 


0.038(0.047) 


0.040(0.156) 


0.094 


0.055 


0.098 



H/3, Mg 2 ,<Fe> 



WITHOUT ENHANCEMENT 

None 20 -0.022(0.057) - - 0.035(0.028) 

All except 2 0.041(0.157) - - 0.039(0.128) 

H/3, [MgFe] 



Table 12. Correlations with log(Age). The number of galaxies 
fitted (N) and best fit slopes and intercepts (errors in brackets) 
are given. Also shown are the (unweighted) correlation coefficient 
(r) and \ 2 values. 





N 


Slope 


Intercept 


r 


x 2 


E/SOs 


[Fe/H] 


17 


-0.701(0.167) 


0.771(0.139) 


-0.81 


71 


[E/Fe] 


17 


0.346(0.116) 


-0.030(0.099) 


0.78 


65 


[Z/H] 


17 


-0.551(0.243) 


0.893(0.202) 


-0.66 


70 


Bulges 


[Fe/H] 


11 


1.658(1.604) 


-0.582(0.814) 


-0.01 


93 


[E/Fe] 


11 


0.364(0.253) 


-0.039(0.11) 


0.34 


45 


[Z/H] 


11 


1.538(1.410) 


-0.387(0.702) 


0.06 


55 



ity dispersion and only marginal (approximately constant) 
enhancement with respect to solar, i.e. [E/Fe] ~ +0.15. 
Thus our results suggest that while the principle parameter 
controlling abundances and abundance ratios in early-type 
galaxies is age, in bulges it is central velocity dispersion. 

We have compared the trends in our early-type galaxy 
sample with the results of Trager et al. (2000b) and 
Kuntschner et al. (2001). Both authors find an anti- 
correlation between age and metallicity that is approxi- 
mately aligned with the age/metallicity degeneracy. These 
correlations are in good agreement with the anti-correlation 
found in our early-type sample (Fig. |l2|). Trager et al. 
(2000b) also found an offset of the age— metallicity relation 
with lower velocity dispersion. Future spectral observations 
of low velocity dispersion ellipticals with high signal-to-noise 
are needed to further investigate these trends, incorporating 
large numbers of indices. 



Table 13. Correlations with log(<ro). The number of galaxies 
fitted (N) and best fit slopes and intercepts (errors in brackets) 
are given. Also shown are the (unweighted) correlation coefficient 
(r) and \ 2 values. 





N 


Slope 


Intercept 


r 


x 2 


E/SOs 


Log(Age) 


17 


1.655(0.690) 


-3.079(1.642) 


0.74 


110 


[Fe/H] 


17 


-0.727(0.881) 


1.906(2.089) 


-0.55 


159 


[E/Fe] 


17 


0.422(0.461) 


-0.735(1.091) 


0.39 


117 


[Z/H] 


17 


-0.553(0.688) 


1.736(1.631) 


-0.53 


31 


Bulges 


Log(Age) 


11 


0.811(1.053) 


-1.265(2.277) 


0.54 


38 


[Fe/H] 


11 


1.043(0.495) 


-2.023(1.093) 


0.80 


49 


[E/Fe] 


11 


0.246(0.423) 


-0.401(0.917) 


0.04 


25 


[Z/H] 


11 


1.177(0.338) 


-2.169(0.750) 


0.85 


15 



4.7 Results from red indices 

Estimation of age and metallicity from the NIR data is dif- 
ficult, as the only age sensitive indices are in the blue wave- 
length range. We must therefore combine indices from pop- 
ulations that may differ, e.g. H/3-CaT in Fig. [Io| The NIR 
indices were also not modelled by TB95, so abundance ratio 
effects can not be estimated. We have nonetheless derived 
log(Age) and [Fe/H] estimates from the calcium triplet in- 
dices and H/3. The Mgl index was excluded as, contrary to 
expectation, the positions of values of this index with respect 
to the SSP grids imply lower metallicities than do values of 
the CaT index. Many of the early-type galaxy data points 
also imply ages greater than 17 Gyr. We therefore assume 
that this index is affected by some, as yet unidentified, cal- 
ibration error (either in our reductions or in the SSP esti- 
mates) . Values of log(Age) derived from the calcium triplet 
and H/3 index are higher than those given in Table [ic| by ~ 
0.15 while [Fe/H] estimates are lower by ~0.2 dex. While we 
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observe that these estimates are, qualitatively, in line with 
our results from blue indices, we make no further attempt 
at interpretation of these NIR data. 



Figure available from 



littp: / /www. star. uclan.ac.uk/~rnp/rcsearch. htm 



5 COMPARISON WITH COMPOSITE 
MODELS 

In this section we generate models with composite SFHs and 
compare the predictions with our observations in an effort 
to understand galaxy histories. The composite model code 
was first described in Sansom & Proctor (1998), where we as- 
sumed solar abundance ratios. Here we extend the models to 
incorporate non-solar abundance ratios. Fourteen elements 
(listed in Table 0) are followed self-consistently. These cover 
most of the heavy element mass loss from SNII, SNIa and 
intermediate mass stars. Our models can allow for inflow 
of gas (enriched to the current level of the ISM or of pri- 
mordial composition) into a single zone. The lowest metal- 
licity modelled by W94 is for Z = Z /1OO. Therefore we 
start the models with 10 6 M© of gas containing this small 
amount of metals (Z = 2 x 10 -4 by mass fraction), assum- 
ing solar abundance ratios within this initial metallicity. Low 
metallicity stars in our Galaxy have an increasing excess of 
a-element abundances (as described by Equation ^). We al- 
low for this in the composite models via the denominator 
in the exponent of Equation fj] which is varied by up to a 
factor of 3 for a-elements in the SSP stars. Allowance for 
a-enhanced SSPs at low metallicity did not produce large 
effects in the predicted line-strengths in galaxies, which are 
dominated by higher metallicity stars. We use V99 SSPs and 
the TOOa method to allow for non-solar abundance ratios in 
the composite models, calling on predictions of line-strength 
changes modelled by TB95. The rate of SNIa is a parame- 
ter in our models. In the current models we use a rate of 
3.8 x 1CT 5 SNIa Gyr" 1 Mg 1 . This is approximately that 
inferred in our Galaxy (Timmes, Woosley & Weaver 1995) 
which has an uncertainty of about a factor of two. Larger 
SNIa rates will produce stronger iron sensitive lines. We as- 
sume a Schmidt law with index of one for the star forma- 
tion rate (SFR=Cxgas density), and a Salpeter initial mass 
function (IMF). 



Figure 13. Luminosity-weighted log(Age), [Fe/H], [E/Fe] and 
[Z/H] are plotted against the log of central velocity dispersion. 
Symbols as Fig. 12. 



5.1 Primordial Collapse 

For the primordial model we started with C=4.0 Gyr" 1 , and 
a rapid, enriched inflow rate of 10 7 Mq Gyr -1 , going down 
to a more steady rate of C=0.2 Gyr -1 and zero inflow after 
0.4 Gyr. Star formation is followed up to 1.5 Gyr ago. In a 
previous paper (Proctor et al. 2000) we showed that (assum- 
ing solar abundance ratios) rapid collapse and star formation 
in a primordial gas cloud does not produce strong enough 
metal absorption lines when compared with observations of 
early-type galaxies and spiral bulges. We confirm this re- 
sult here with our non-solar ratio models. This is shown in 
Figs |^ and ^ where the thickest, short line indicates our 
predictions for such a primordial collapse model, for times 
ranging from 10 to 17 Gyr after the start of star formation. 
These primordial models include a higher rate of 7.6xl0~ 5 
SNIa Gyr -1 Mg 1 . The predicted metal line-strengths are 
too low to account for the observed line-strengths. Thus the 
conventional picture of spheroid formation through rapid, 
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early collapse, followed by passive evolution, is ruled out by 
the observed strong lines. A similar conclusion was found 
for nearby spheroids by Worthey, Dorman & Jones (1996) 
and is analogous to the well known G-dwarf problem in our 
own Galaxy, where there are insufficient low metallicity stars 
compared to predictions of closed box models with stars gen- 
erated with a Salpeter IMF. 



5.2 Models with extended inflow 

In Proctor et al. (2000) we found that observations of spi- 
ral bulges could be explained with extended inflow models, 
with gas inflow over several Gyr enriched to the current 
level of the ISM. This assumed solar abundance ratios. Al- 
lowing for non-solar ratios we find that such models tend 
to under-produce Fe sensitive features. This is because the 
early feedback from SNII is extremely Mg rich compared to 
Fe (several 10s of times solar ratios - see the SNII models of 
Woosley and Weaver 1995). To produce models which can 
simultaneously explain Fe and Mg sensitive spectral features 
a delayed burst of star formation seems to be needed to al- 
low the ISM to first become enriched with Fe peak products 
from SNIa. Such models are preliminarily explored in the 
next section. 



5.3 Merger Models 

With a delayed burst of star formation we begin to be able 
to produce models which can simultaneously explain the 
strengths of several spectral features. A full exploration of 
SFH parameter space is beyond the scope of this paper and 
will be the subject of future work. However, in Figs ^| and 
H we illustrate predictions of line-strengths for merger mod- 
els with an associated burst of star formation several Gyr 
after the start of the SFH (delayed burst models - medium 
thick short line). Present day predictions (at 17 Gyr) for 
models with a burst ranging from 3 Gyr to 13 Gyr delay 
are shown. The parameters used to describe the composite 
merger model shown here are C=4.0 Gyr -1 initially, with 



( |1996| ) used a semi-analytic model of galaxy formation in 



both field and cluster environments to make testable pre- 
dictions for the ages of various Hubble types in differing 
environments. In these models early-type galaxies form by 
the merger of two, roughly equal mass, progenitor galax- 
ies, while spiral galaxies form by accretion o f a d isc onto 
a pre- formed elliptical galaxy. Kauffmann (1996) found 



luminosity-weighted ages for early types between 5 an d 12.5 
Gyr, in good agreement with our findings. Kauffmann (199£ ) 
also found cluster ellipticals to be ~ 4 Gyr older than ellip- 
ticals in low density environments. This is again consistent 
with our findings, as the five Virgo cluster ellipticals are 
amongst the oldest in our sample. The positive correlation 
between age and velocity dispersion in elliptical galaxies is 
qualitatively consis tent with the relationship proposed by 
Forbes & Ponman (199£) for a versus merger redshift. Con- 
sequently, while both primordial collapse and extended in- 
flow models fail to reproduce the main features of our data 
without recourse to additional physics (a biased IMF or Pop- 
ulation III stars), the hierarchical model of early- type galaxy 
formation by merger agrees well. 



For the bulges of spirals, Kauffmann (1996) predicts 
ages significantly lower than those in ellipticals , in a gree- 
ment with our findings. However, Kauffmann (1996) also 
predicts a correlation between bulge ages and the luminos- 
ity of their discs. Inspecting the distribution of Hubble types 
(indicated by solid symbol sizes in Figs |l| and |l|) we find 
no evidence of a trend in age with spiral Hubble type. There- 
f ore, t his prediction is at odds with our findings. Kauffmann 
( 1996| ) does point out, however, that the correlation may 



be hidden if there is significant inflow of gas from the disc, 
perhaps due to the formation of bars. Many of the bulges in 
our sample show strong evidence for kinematic sub-structure 
and on-going star formation (i.e. emission). Our bulges also 
possess [Mg/Fe] values too low to have been formed in a 
single primordial burst (see Fig. ||). These observations, and 
the absence of a correlation between age and Hubble type, 
support the idea that disc inflow must play an important 
role in the star formation histories in bulges. 



10 7 Mq Gyr -1 during the star- burst. This rapid star-burst 
lasts for 0.4 Gyr after which the inflow is set to zero, so the 
remaining gas is rapidly used up in star formation. We see 
that these examples reach the regions populated by early- 
type galaxies, for most Fe and Mg sensitive line-str engths . 
An exception is Fe5406 (but see Fig. ^ and Section 3.2.1). 
Interestingly, we find that the models tend to under-predict 
carbon sensitive features (<CN>, C24668). This may indi- 
cate that all the sources of carbon enrichment have not been 
accounted for in our models. Indeed, carbon enhancements 
in some very low metallicity stars in our own Galaxy are 
hard to explain (Norris, Ryan & Beers 1997) and dredge- 
up models for the contributions of carbon from intermedi- 
ate mass stars are uncertain. Stars produced in the delayed 
burst have lower overall [E/Fe] than the earlier star forma- 
tion, since SNIa have had time to accumulate Fe in the ISM. 
Thus these model predictions tend to support the idea of el- 
lipticals forming by mergers of galaxies, with enriched gas 
inflow and enhanced star formation during the merger. 

There are few predictions of detailed galaxy proper- 
ties from hierarchical merger models. However, Kauffmann 



6 CONCLUSIONS 

We have derived luminosity- weighted log(Age), iron abun- 
dance, abundance ratios and metallicity (log(Age), [Fe/H], 
[E/Fe] and [Z/H] respectively) in the centres of 32 galaxies, 
ranging in Hubble type from E to Sbc. We used 20 indices 
for which V99 modelled SSPs, plus their sensitivities to indi- 
vidual elements as tabulated in TB95, to model the effects of 
non-solar abundance ratios in these galaxies. We find that ig- 
noring such enhancements leads to reasonably accurate age 
and metallicity estimates, but that using many fewer indices 
(e.g. 3) leads to larger errors in derived values (see Table |ll|) . 
By using many indices and modelling abundance ratios we 
are able to probe correlations between derived parameters 
much more accurately. 

Our sample of early-type galaxies spans a wide range 
of ages (2 to 13 Gyr) with the 5 Virgo cluster ellipti- 
cals amongst the oldest. These E/S0 galaxies show corre- 
lations of both velocity dispersion (a) and [E/Fe] with age, 
while [Fe/H] and [Z/H] both show anti-correlations with age. 
[Fe/H], [Z/H] and [E/Fe] show no significant correlation with 
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a over the small range in a covered by our E/SO galaxies. 
These results are at odds with the predictions of primor- 
dial collapse models, which predict uniformly-old, early-type 
galaxies and increasing [Z/H] with a. However, the correla- 
tions suggest that the main parameter controlling the metal 
content of our bright, early-type galaxies is age. Our results 
are consistent with the predictions of hierarchical merger 
models of galaxy formation (e.g. Kauffmann 1996) which 
predict the observed E, SO (and bulge) age sequence. The 
observed correlation between age and velocity dispersion is 
in qualitative agree ment with the relationship proposed by 
Forbes & Ponman (|L999|) for merger remnants. The strong 
correlation of [E/Fe] with log(Age) and anti-correlation of 
[Fe/H] with log(Age) are qualitatively reproduced by our 
merger models of galaxy formation (Section and Fig. |^). 
In our models these trends are the result of the shortening 
interval (in which SNIa can produce Fe) between the com- 
mencement of star formation and the final merger event. 
Thus we find several observational results which agree with 
the predictions of merger models for early-type galaxies. 

The correlations outlined above in turn call for careful 
interpretation of Lick index versus a correlations. For in- 
stance, high Mg2 , high velocity dispersion, early-type galax- 
ies are amongst the most metal-poor in our sample. Any 
correlation of Mg2 with a, may in fact, reflect an increas- 
ing age and enhancement ([E/Fe]) with velocity dispersion 
(indicating mass) not, as is usually assumed, a metallicity- 
mass relation (s ee also related comments in the conclusion 
of Trager et al. ( |2000b| ) ) . 

The anti-correlations of [Fe/H] and of [Z/H] with 
log(Age) could also have implications for the interpreta- 
tion of colour-magnitude diagrams of these objects. This 
can be illustrated by the extremely narrow range of opti- 
cal colours (A(U-V) ~ 0.08 for SSPs ranging in age from 3 
to 17 Gyr) predicted for SSPs that follow the age-[Fe/H] 
anti-correlation shown by our early-type galaxy sample. 
Colour-magnitude correlations have also been interpreted as 
metallicity-mass relations on the basis of monolithic collapse 
models for the formation of elliptical galaxies (e.g. Kodama 
et al. 1998). Our results from line-strengths suggest that 
merger models must be considered before this interpreta- 
tion of the observed correlations can be relied on. 

We detect significant differences between early-type 
galaxies and spiral bulges. Es, SOs and bulges in our sam- 
ple form a continuous overlapping sequence of decreas- 
ing luminosity- weighted age, with bulges typically 2 Gyr 
younger than SOs, and 5 Gyr younger than Es. T his is again 
in line with the predictions of Kauffmann (1996), for galax- 



ies in low density environments. We find no significant cor- 
relations with age in bulges. However, correlations of [Fe/H] 
and [Z/H] with a are strong. Thus the main paramete r con - 
trolling the metal content of bulges is a. Kauffmann (199C) 
predicts a correlation between bulge-to-disc ratio and age in 
spiral galaxies. We de tect n o such correlation in our data. 
However, Kauffmann ( 1996[ ) also points out that if signifi- 
cant inflow from the disc occurs, after a merger, this cor- 
relation may be lost. The relatively low values of [E/Fe] 
(~ +0.15) in bulges and the strong correlations of metals 
with a are consistent with such a model of bulge forma- 
tion. The correlations of [Fe/H] and [Z/H] with a in spiral 
bulges means that we interpret the correlations of Lick in- 
dices with a found in this work (Figs ^ and ^ and Table [5]) 



as a metallicity-mass relation, in contrast to our finding for 
early- type galaxies. 

In conclusion, we have shown that primordial collapse 
models of galaxy formation are unable to reproduce the line- 
strengths observed in the spheroids of galaxies, while merger 
models can. Derived ages and correlations between derived 
parameters differ significantly between early and late-type 
galaxies, suggesting that, at least at some point in their 
evolution, the star formation histories in these objects must 
have differed significantly. We therefore contend that the 
similarities in morphology and photometric properties in 
these objects, noted in the introduction, are the result of 
the various degeneracies at work rather than indicating sim- 
ilar formation processes. 
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APPENDIX A: VELOCITY DISPERSION 
CORRECTIONS 

Al Characterisation of stellar indices with 
spectral broadening 

For accurate calibration, indices in galaxy spectra require 
correction to account for the effects of internal stellar veloc- 
ity dispersion. Consequently, the indices of a sub-sample of 
Lick calibration stars were measured after convolving their 
spectra with Gaussians of a range of widths. 15 stars (spec- 
tral types G8 to K3), which best matched the spectral energy 
distributions of the galaxies, were used for the blue indices. 
In the red, all 24 Lick calibration stars were found to match 
galaxy spectra well. For each index a correction factor was 
calculated at each value of broadening. For most line indices, 
the correction factor was calculated as the ratio (Ci) of the 
index value at the Lick resolution to that at each value of 
broadening, i.e: 



Ci — Ih/lMeas 



(Al) 



Where I_l is the index value at the calibration resolution 
and Lifeas is the index value in the broadened spectrum. 

For molecular band indices (CN1, CN2, Mgi and Mg 2 ) 
and indices with ranges spanning zero (H5 and H7 indices), 
correction factors were calculated as the difference between 
the measured index and that at the calibration resolution at 
each value of broadening, i.e: 



Ci — II — Ihleas 



(A2) 



The appropriate correction factor was calculated for 
each index, at each value of broadening, by averaging values 
from the stellar sub-samples. A polynomial fit of order 3 was 
found such that: 

Ci — xo + x\oc + X20~c + £300- (A3) 

Where oc is the width of the Gaussian (in km s _1 ) 
convolved with the stellar spectrum at ol- Behaviours of 
the indices of the stellar sample with broadening above the 

as lines. These are in 



Al 



Lick resolution are shown in Fig. 
good agreement with the behaviours of indices measured 
by previous authors (e.g. Kuntschner 2000). The nature of 
the applied velocity dispersion corrections, i.e. whether Ci 
should be multiplied by the raw index value or added to 
the raw value are indicated in the plots and summarised 
in column 2 of Table |^]. This table also gives a typical 
uncertainty in the correction factor for each index, estimated 
as the RMS scatter in the stellar values about the mean 
Ci at o (7=200 kms . This corresponds to the appropriate 
value for a high velocity dispersion elliptical galaxy. When 
correction factors are used the scatter in the stellar data 
at the value of broadening appropriate to the galaxy was 
included in the statistical errors. 



differences in the behaviour of galaxies compared to the stel- 
lar data were observed. These indices (most noticeably H<5, 
Ca4227 and Ca4455) are among those affected by the poor 
removal of the dichroic response. The scatter in the stellar 
dat a has been allowed for in our errors as described in Sec- 
tion 2.4.2. However, to test the possible impact of differences 



in behaviour between our samples of stars and galaxies, two 
velocity dispersion correction techniques were tested. 

Technique 1. 

The first technique (which is generally used by other au- 
thors) is to broaden all galaxy spectra by a Gaussian of a 
width ob, given by Equation [j], prior to measurement of 
indices. The resultant spectrum of a galaxy with velocity 
dispersion ay, then has a total broadening (og) given by: 



2 

Oq 



2 
°7 



2 



2 . 2 

Oy + o B 
2 

f o v 



(A4) 
(A5) 



As the correction factors are dependent on the excess 
broadening of the galaxy with respect to ox the value oc 
— oy is entered into the polynomials in Table Al to obtain 
the correction factor Ci. 

Technique 2. 

In the second technique, if the total broadening of the galaxy 
(oj+oy) 1 / 2 is less than the target calibration resolution 
(ol), the galaxy spectrum is broadened up to the Lick res- 
olution by convolution with a Gaussian of width ob given 
by: 



2 

o B 



2 

: o L 



2 
Oj 



2 

o v 



(A6) 



Consequently, as the spectrum has been broadened to 
the calibration resolution, the final index values can be mea- 
sured directly from the broadened galaxy spectrum with no 
need for correction. 



If, on the other hand, the total galaxy broadening is 
greater than ol, then the spectra are left un-broadened and 
oc=(oj + o v — Ol) 1 '' 2 is substituted into the polynomials 
to calculate the appropriate value of Ci. 



It should be noted that, for some indices (e.g. H5), the 
second technique involves no use of the polynomials as all 
galaxies in our sample have o\ > oj + o v . The two tech- 
niques agree well over our whole velocity dispersion range. 
This is true for all indices measured. Comparisons of the re- 
sults from the two techniques can be found in Proctor (PhD 
thesis - in preparation). In this paper the second technique 
is used since it involves smaller corrections, thus minimising 
the effects of systematic differences in the behaviour of the 
stellar and galaxy spectra with broadening. 



A2 Application to galaxies 

To gain confidence in these polynomials, galaxy spectra were 
also broadened by convolution with Gaussians of a range of 
widths. The behaviour of the galaxies was then compared to 
that of the stellar sub-sample (Fig. Al). For most indices the 
galaxy data lie within the scatter of the stellar data. How- 
ever, for a few indices, large scatter and/or small systematic 
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Figure available from http://www.star.uclan.ac.uk/ rnp /research. htm 



Figure Al. Behaviour of indices with broadening. Stellar correlations are shown as lines. Galaxy data are shown as points. Only galaxies 
with total observed broadening below ax are plotted. 
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Table Al. Polynomial coefficients relating the velocity dispersion correction factor (C;; see text) to the correction velocity (cc). 



Index 


(A)dditive 
or 

(M)ultiplicative 


Polynomial Coefficients (Cj 

XQ Xl X2 

(x to- 3 ) (x icr 6 ) 


versus in km s 1 ) 
(x 10- 9 ) 


Uncertainty 
in correction 
at cr c =200 km s" 1 


H<5 A 


A 


0.0 


-0.058 


-5.195 


3.894 


0.029A 


H<5 F 


A 


0.0 


-0.005 


0.704 


-0.077 


0.019A 


CNi 


A 


U.U 


U.OUz 


O.Uy ( 


-0.081 


O.OOlmag 


CN 2 


A 


0.0 


0.003 


0.225 


-0.201 


O.OOlmag 


Ca4227 


M 


1.0 


0.231 


1.246 


13.400 


3.3% 


G4300 


M 


1.0 


0.021 


0.480 


0.133 


0.4% 




A 


0.0 


0.006 


1.181 


-8.883 


0.053A 


H7F 


A 


0.0 


-0.079 


-2.751 


4.096 


0.054A 


Fe4383 


M 


1.0 


0.037 


2.334 


-0.811 


1.7% 


Ca4455 


M 


1.0 


0.086 


3.941 


1.424 


4.1% 


Fe4531 


M 


1.0 


0.033 


1.608 


0.628 


1.4% 


C 2 4668 


M 


1.0 


-0.001 


0.851 


0.284 


0.5% 


H/3 


M 


1.0 


0.031 


0.213 


0.384 


1.0% 


Fe5015 


M 


1.0 


0.040 


3.025 


-2.443 


1.1% 


Mgi 


A 


0.0 


0.002 


0.067 


-0.073 


O.OOlmag 


ivig2 


j\ 


0.0 


0.004 


0.035 


-0.008 


O HOI mno- 


Mgb 


M 


1.0 


-0.003 


2.129 


0.307 


0.4% 


Fe5270 


M 


1.0 


0.024 


3.027 


-2.352 


0.5% 


Fe5335 


M 


1.0 


0.059 


4.748 


1.522 


0.6% 


Fe5406 


M 


1.0 


0.050 


3.927 


3.565 


1.2% 


Cal 


M 


1.0 


-0.053 


1.892 


2.526 


0.1% 


Ca2 


M 


1.0 


0.004 


1.248 


2.339 


0.1% 


Ca3 


M 


1.0 


0.028 


1.592 


0.909 


0.1% 


Mgi 


M 


1.0 


-0.082 


11.286 


-1.443 


0.1% 
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